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, Abstract. We construct differential invariants that vanish if and only if the geodesic flow 

' of a 2-dimcnsional metric admits an integral of 3rd degree in momenta with a given Birkhoff- 

, Kolokoltsov 3-codifferential. 
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1. Introduction 



1.1. Definitions and results. Let Shea, surface (i.e., 2- dimensional real manifold) equipped 
with a Riemannian metric g given in local coordinates hj g = • gijdx^dxK Since the metric 
g allows us to identify the tangent and cotangent bundles of S, we have a scalar product 
Q \ and a norm on every cotangent plane. The geodesic flow of the metric g is the Hamiltonian 

■ system with the Hamiltonian H := = ^g^-'PiPj, where p = (pi,P2) are the momenta and 
J3 ' M is the norm induced by g. 

We say that a function F : T*^ M is an integral of the geodesic flow of g cubic in 
momenta (shortly: cubic integral for the metric g), if 

J> ' (1) in the local coordinates x := x^,y := x"^ on the surface and the corresponding momenta 
^ ■ Px,Py, -F is a homogeneous polynomial in the momenta of degree 3: 

^ : F{x,y;p^,py) = ao{x,y)pl + ai{x,y)plpy + a2{x,y)p^pl + a3{x,y)pl, (1.1) 

0> ' (2) F is an integral of the geodesic flow of g, i.e., {H, F} = 0, where {■, ■} is the canonical 
^ ■ Poisson bracket on T*S. 

^ I Cubic integrals allow to construct different geometric structures on the surface; we shall 
use the so called holomorphic BirkhofF-Kolokoltsov 3-codifferential . We consider the complex 
^ ' structure on S given by g: a local complex coordinate z = x + \y is determined by the 

■ property that the metric g has the isothermic form 

g = A(x, y){dx'^ + dy'^) = \(z, z)dzdz. (1.2) 

Then, the Birkhoff-Kolokoltsov 3-codifferential (corresponding to the cubic integral (11.11) ) is 
a complex- valued symmetric (3, 0)-tensor given by 

A{z) := (ao(x, y) - a2(x, y)) + i (ai(x, y) - a^{x, 2/)) ) f ® £ ® £• (1-3) 

As it was shown by Birkhoff [Bi], the coefficient (ao — 02) + i (ai — 03) is a holomorphic 
function of the complex coordinate z = x + \y. Kolokoltsov observed in [Kol] that under a 
holomorphic ( ^ conformal) coordinate change this coefficient changes as the corresponding 
coefficient of a complex (3, 0)-tensor. Below in § 12.21 we shall prove both properties of A and 
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give several useful formulas for A in arbitrary (not necessarily isothermic) coordinates, such 
as how to compute A from F and how to check the holomorphicity of A, see Lemma 15. 1[ 
The goal of this article is to give an answer to the following question: 

Given a metric g and a holomorphic 3-codifferential A, how to decide whether there 
exists a cubic integral for the geodesic flow of g whose Birkhoff-Kolokoltsov 3-codif- 
ferential coincides with A ? 

If such a cubic integral F exists, we say that A is compatible witli g. 

A complete algorithmic answer is given in § 11.31 using Theorems II. ![ 11.21 and also Proposi- 
tions 13.11 14. 2[ Moreover, in the most interesting case. Theorems 11.11 11.21 give an explicit 
formula for the cubic integral. 

1.2. Main theorems. In the formulas below we use Einstein's summation convention, ie., 
we sum over repeating upper and lower indices. The notation a-jk and so on means covariant 
derivation(s) of a function (or a tensor) a, and a^*-''^^ denotes the symmetrisation in indices 
ij,k. 

Let X = x^,y = x"^ he arbitrary coordinates on the surface 5*. Given a metric g = 
gijdx^dx^ on S, let Rljj^ be its Riemannian curvature, and set X{x,y) := ^ydet{gij). Since 
our considerations are local, we equip 5* with the orientation given by the metric volume 
form Ug = Xdx A dy. 

We denote by {f,h}g := ^ (^f^g - the Poisson bracket of functions f,h on S 

with respect to Ug. Further, the tensor J^j := g'^^Ukj is the operator of rotation by 90°, it 
defines the complex structure on S compatible with the metric g and the orientation, i.e., 
corresponds to the multiplication by i in the complex coordinate z = x -\- '\y. 

Define the following smooth functions y^o, . . . , 953, which are invariant algebraic expressions 
of the components of g and its derivatives: 

• The Gauss curvature, which is the half of scalar curvature: 

:= i? := (1.4) 

• Half-square of the gradient of ip^ = R, half-square of the gradient of the result: 

:= i|Vv.oP = l^^-'S^lf ; ^3 := ||V^,p = Ig-'^'^ (1.5) 

• The Poisson bracket of (po, ipi with respect to Ug-. 

^,:={^o,V^i}.:=i(^^-^t^)- (1-6) 

Further, let A be a complex 3-codifferential on 5" given in some (not necessarily isothermic) 
coordinates (x^,x^) as symmetric (3, 0)-tensor A^^'' with complex components. Define A : = 
^{A) to be its real part, or A^^^ := ^{A^^'' + A^^'') on the level of components where a 
means complex conjugation. Then A = [A^^'^) is a (usual) symmetric (3, 0)-tensor with real 
components. The formula for the imaginary part '^{A^^'') is given in Section [5l 

In addition to (/jq, • • • , V^s, we define the functions Dq, . . . , D3, ^2, ^3, 1-form /C = ]Cidx\ 
and symmetric (3, 0)-tensors B^^^, F^^^. They all are certain invariant algebraic expressions 
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in g, A, and their derivatives. 

A :=5^^(Do),(^o);, - 4i*^^, ■ (ifo),, ■ {ifo),k; (1.8) 
D2:={Do,Vi}g + Wo,Di}g- (1.9) 
Ds ■.=g'\Di),,{^^),, -AA'^K^, ■ {^^).^ ■ (<^i);fc (I-IO) 

All these functions have a clear geometric sense: 

• Dq is the triple divergence of 4 ■ A; 

• g^^ {Do)-i{ipQ).j is the scalar product (V-Dq, Vy^o) of gradients, 

• A'^^^.i{ipo)-j{(po)-k is the evaluation of the quadratic form corresponding to symmetric 
(2,0)-tensor {d'wAy^ := A'^^'-i on the 1-form dipo = dR; 

• We obtain (p^ from (fi by the same formula as (fi from ifQ, similarly, we obtain from 
Di and (pi by the same formula as Di from Dq and </?o- 

Finally, define 

^2 := T^etl (v^i);.. {Vi);y Di\ := ^ detl {^i),, {^i),y Di \ , (1.11) 





/Ci := — detf^'^^J'^ B'^^ := g'^'^g^^^J'^ilCm p'^^ := A'^^ + B'^^ . (1.12) 



Theorem 1.1. /n the above notation, assume that the 3-codifferential A is holomorphic 
and that >p>2 in non-vanishing on S . Then there exists a cubic integral F with the Birkhoff- 
Kolokoltsov 3-codifferential A if and only if g and A satisfy the equations Q2 = 0, Q3 = 0. 

Moreover, such F is unique and in local coordinates {x^,x'^) with the dual momenta (pi,P2) 
the integral F is given by 

F{x^,x'^;pi,P2) := F'^''{x^ , x'^)piPjPj 
where the symmetric {3, 0) -tensor F^^'^ is computed using the formulas (11.41) - fll.l2p . 

Remarks. 1. Our definitions and formulas involve the orientation on the surface S, however, 
the final result is, of course, independent of the orientation. More precisely, after reversing 
the orientation every real expression (formula) either remains unchanged, or inverts the sign. 

2. After appropriate "cosmetic" changes. Theorems 11.11 [L2] remain valid also in the pseudo- 
Riemannian case, see Proposition 16. 2[ 

Let us now consider the degenerate case when ip2 = {R, ||V-Rp} vanishes identically. In 
this case we set ip^ := ipo, := Dq 

y^l := AR DI := AD^ - 23ft((A,i?;,);,) . (1.13) 

and define the ^-versions ip2,f3,G2^D2,D^, Q^, ICl, IC2.-. of the previous expressions by the 
replacing (/^o, ^0, V'l, - by y?*, D*,ipl,... in the formulas (USD, ([L9]), ffLTOD . ffTTTD . 
above, see also Section HI Further, let V be that of two expressions Vi := detQ^^"!'""' j 

(non-*-versions!) for which (v7o);x» is non-vanishing (any of two if both {ipo)-x, {'^o);y are 
non-zero.) 
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Theorem 1.2. Let g be a metric and A a holomorphic 3-codifferential. Assume that the 
function (p2 = {R, ||V-R|^}g vanishes identically and the function = {R, AgR}g is non- 
vanishing. 

Then g admits a cubic integral F with the Birkhoff-Kolokoltsov 3-codifferential A if and 
only if g satisfies the PDEs ^25^3; ^'^'^ Moreover, such F is unique and in local coordi- 
nates with the dual momenta (^1,^2) ihe integral F is given by 

F{x\x^;pi,P2) := F'^''{x\ x^)piPjPj 

where the symmetric {3,0)-tensor F^^^ is computed by the *— version of the formula (11.121) . 
namely, 

K* := ^det(^|pj'^ B^^^ := g^'^ g^^^ r\]C^ F'^^ := i*^' + . (1.14) 

Remark. The equations ^3 (always) necessary conditions for compatibility of A 
with (7, but not sufficient in general. Nevertheless, if A is compatible with g and is non- 
vanishing, then the formulas (11.141) are valid and give the correct value of the cubic integral 
F, even in the case when is also non-zero. 

The only two cases which are not covered by Theorems 11.11 11.21 are the case ip^ = const 
(i.e., the metric has constant curvature) and the case ip2 = 0, V92 = 0. These cases are easier, 
and will be solved in Propositions 13.11 14. 2[ 

I. 3. Algorithm for checking the existence of cubic integrals. Let g he a, metric on 
a surface S and A a holomorphic 3-codifferential. Theorems 11.11 II. 2[ classical results of 
Darboux and Eisenhart which we recall in §1.5^ and also Propositions 13.11 14.21 give us the 
following algorithm to decide whether there exists a cubic integral F with the given Birkhoff- 
Kolokoltsov 3-codifferential A. Moreover, in the most interesting cases covered by Theorems 

II. H 11.21 we obtain an explicit formula for F as algebraic expression in (7, A and their deriva- 
tives. In other cases covered by Propositions 13. 1| 14. 2| the formula for F requires a solution of 
linear systems of ODEs. Corresponding calculations can be easily realized using computing 
software such as Maple® or Mathematica®. 
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Let us comment some boxes in the flowchart of the algorithm. Recall that A is compatible 
with (yf if A is the Birkhoff-Kolokoltsov 3-codifferentials associated to some cubic integral F 
of the metric g. 

(bl) If R is constant, then every cubic integral is of the form F = o^^^LiLjLk with constants 
Q^jk ^ where Lj, z = 1, 2, 3 are three linear (polynomial of degree 1) invariants which 
are linearly independent. Thus the cubic integrals form a linear space of dimension 10, 
see |Krj . On the other hand, the space of holomorphic 3-codifferentials A has infinite 
dimension. Proposition 13.11 shows that the equation Z^o = is necessary and sufficient 
condition for compatibility of A with g and, if it is fulfilled, the space of cubic integrals 
depends on 3 real parameters. More precisely, the a generic cubic integral has the form 
3f?(A) + Upi ■ H where h = (6*) is a vector field with given values b{P) and rot b{P) at a 
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given point P E S. The construction of the vector field b = (6*) reduces to solution of 
certain linear systems of ODEs. 

(b2) It is sufficient to calculate only one of the expressions Vx,Vy, namely one such the 
corresponding is not zero, see Theorem 11.21 

(b3) By a classical result of Darboux and Eisenhart (see discussion in §1.5p the vanishing 
of both and (/?2 is equivalent to the local existence of a Killing vector field which 
we denote by C,- The latter is equivalent to the existence of a linear integral which we 
denote by L. Since Lie derivative of the curvature R along ^ must vanish, we have that 
^ is proportional to the Hamiltonian vector field of i?, ^ = f{x,y) ■ Xr. Such f{x,y) 
is unique up to a constant factor, and can be found from the PDE-system ^^^^g = 
on the unknown function /. Note that the PDE-system J^j^^g = is Frobenius and, 
in view of conditions 9?2 = and = 0, is in involution, i.e., is essentially a system of 



(b4) The expressions and the part of the algorithm starting from the box (b3) is 

treated in § I01 

In the pseudo-Riemannian this algorithm basically works, however there exists an impor- 



consequence we obtain a new special class of pseudo-Riemannian metrics whose counterpart 
does not exists in the Riemannian case. They are not metrics of constant curvature, they 
have ip2 = = ip2, however, as we prove in Proposition 16.31 the metrics from this class admit 
no cubic integrals, and therefore no Killing vector fields. We refer to Section [6] for further 
details and description of an additional step in the algorithm. 

We emphasise here that the algorithm works in any (not necessary isothermal) coordinate 
system. In particular. Lemma [5.11 gives a way how to verify whether the data {g,A) have 
sense, i.e., whether real (3,0)-tensor A = ^{A) corresponds to the real part of a holomorphic 
3-codifferential. 

1.4. Motivation and history. The importance of polynomial integrals for studying the 
geodesic fiow was recognised long time ago. Indeed, it was Jacobi who had realised that the 
geodesic fiow on the ellipsoid admits an "extra" quadratic integral. This allowed Jacobi to 
integrate the geodesic flow on the ellipsoid. 

Polynomial integrals for the geodesic flow are interesting for physics and for differential ge- 
ometry. The interest from physics is due to the following observation (called the Maupertuis 
principle, or coupling constant transform): for sufficiently large energy levels h the restriction 
of the Hamiltonian system on T*M with the Hamiltonian H := + (where is again 
the kinetic energy corresponding to g, and ^ : S* ^ M is the potential) to the energy surface 
{(x,p) G T*M I H{x,p) = h} has the same unparametrised orbits as the geodesic fiow of 
the metric g^ := {h — V)g. Moreover, if the Hamiltonian system is polynomially-integrable, 
then the geodesic fiow is also polynomially-integrable. More precisely, if a function 



ODE. 



tant difference. Namely, the condition |V-R 



does not imply that R = const. As a 



F{x, y; p^, py) = aopl + aiplpy + OaP^pJ + a^p^ + Cop^ + CiPy 



(1.15) 




Fz{x,y;p^,Py) Fi(x,y;px,py) 



is an integral of the Hamiltonian system with the Hamiltonian + V, then the (homo- 
geneous cubic in momenta) function 

Ffr-= Fs + Hh ■ Fi (1.16) 
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is an integral for the Hamiltonian system with the Hamiltonian ■= 2{h-v) 1^^' 
the geodesic flow of the metric gt '■= {h — V)g. 

Note that the assumption that the integral fll.lSp does not contain the terms Fq and F2 
of degree and 2 in momenta is a natural one: by Whittaker [Whj . if a function of the 
form Fq + Fi + F2 + F3 is an integral for the Hamiltonian system with the Hamiltonian 
H := + V, then Fi + F3 is also an integral. 

Metrics admitting integrals polynomial in velocities are also interesting for geometry, since 
the existence of such integrals implies interesting geometric properties of the metric. More- 
over, study of polynomial integrals helps to solve pure geometrical problems. For example, 
the existence of the integral of degree 1 is equivalent to the existence of a Killing vector field, 
and implies the existence of the coordinate system such that the component of the metric 
does not depend on one of the coordinates [Da] . The existence of an integral of degree 2 is 
equivalent to the existence of another metric having the same geodesic with the initial one 
[Ma-Toli IMa-To2j . The last observation was recently used |Br-Ma-Mal IMa] in the solution 
of two geometric problems explicitly formulated by Sophus Lie in 1882. 

It is not clear whether the existence of the integral of degree 3 has such a clear geo- 
metric meaning, but still, as we mentioned above, it allows one to construct a holomophic 
3-codifferential with interesting properties, and, as it was shown in |Du-Ma-Tb] . a volume- 
preserving vector field on the surface. 

1.5. Previous results in this field. The most classical result in this direction is in folklore 
attributed to Bonnet, and appeared in the books of Darboux [Daj and of Eisenhart [Ei] . 

Given a metric g = {gij{x, y)) on a surface, let us consider the functions ipo, ^i, ^i, f2, f2 
from §1.11 (fQ := R is the Gauss curvature, ipi := ||V-R|^ is half-square of the gradient of R, 
ipl := AgR = Agipo the metric Laplacian of i? = ipo, and ip2 := {v^o, '^i}g, ¥^2 •= {v'o, V^i} are 
the Poisson brackets. 

One can show (see Darboux [Pat §§688, 689] and Eisenhart [Eil PP- 323-325]) that the 
metric admits a nontrivial integral linear in momenta (locally, in a neighbourhood of almost 
every point) if and only if the invariants (f2, ^2 vanish. 

Though Darboux and Eisenhart proved this result for Riemannian case only, their proof 
can be generalised to the pseudo-Riemannian case with the following reformulation of the 
assertion: if the equations (p2 = and ip2 = are satisfied and ||V-R|^ = ipi is non-zero, 
then g admits a nontrivial linear integral. We study the class of pseudo-Riemannian metrics 
with the condition |V-R|^ = in Section [6l and prove that they admit a linear integral iff 
R = const. 

We see that we have an algorithmic way to decide whether a given metric admits a linear 
integral of degree: one need to calculate the invariants ip2, ^2 which are algebraic expressions 
in the components of the metric and their derivatives up to order 5, and compare them with 
zero. 

It was a long standing problem to find a similar algorithmic way to decide whether a 
given metric admits an integral quadratic in momenta. The first attempts are due to Roger 
Liouvill^ [Li] and Koenigs [Koej . They understood that the PDE-system for the coefficients 
of the integral is linear and of finite type, implying that that there must exist an algorithmic 



Roger Liouville was a younger relative of the more famous Joseph Liouville and attended his lectures at 
the Ecole Polytechnique 
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way to decide whether a given metric admits an integral quadratic in momenta. Unfor- 
tunately, the calculation were too complicated to be fulfilled in that time. An interesting 
approach to this problem is due to Sulikovski |Suj : he invented a trick which allowed to 
simplify the calculation. Finally, the problem was solved in recent independent works by 
Kruglikov [Kr] and Bryant et al |Br-Du-Ea] . 

Geodesic flows admitting cubic integrals is a much more difficult subject than those ad- 
mitting linear or quadratic integrals. Indeed, metric admitting linear and quadratic inte- 
grals are completely described: the local description is done by classics jPal §§540-596], see 
|Bo-Ma-Pul[ IBo-Ma-Pu2] for the discussion of the pseudo-Riemannian case, and the global 
(= if the surface is closed) is due to Kolokoltsov |Kolj and Kiyohara [Kilj . see the survey 
|Bo-Ma-Fo] . The metrics admitting linear integrals are parameterised by one function of 
one variable, the metrics admitting quadratic integrals are parameterised by two functions 
of one variable. No such description is known in the cubic case (though from general theory 
of PDE it follows that locally, in the analytic formal category, the metrics admitting cubic 
integrals are parametrised by 3 functions of one variable). 

Though the complete description of metrics admitting cubic integrals is not known yet, 
there exists a lot of local examples (the first local examples appeared already in Darboux 
|Dal §619]). Some of most interesting examples come from physics, see the survey [Hi], or 
are motivated by physics |Hal lYel IVe-Ts] . 

It appears to be difficult to construct metrics admitting cubic integrals on closed surfaces. 
As it was shown by Kozlov [Koz] and Kolokoltsov |Kolj . the surfaces of genus ^ 2 do not 
admit (nontrivial) polynomial integrals; we repeat the argument of Kolokoltsov in the proof 
of Corollary 12.11 Of course, the metrics admitting linear integral also admit cubic integrals, 
since the cube of the linear integral is a cubic integral. Besides these half-trivial examples, 
the list of known examples is very short: there is a series of explicit examples coming 
from Goryachev case of rigid body motion and its generalisation, a series of recently found 
examples due to [Du-Ma] . and proof of the existence of other examples in [Ki2] . [Se] . 

In particular, in all known examples the surface is the sphere (i.e., there is no known 
example of metric on the torus admitting cubic integral and do not admitting linear integral; 
moreover, Bolsinov et al conjectured |Bo-Ko-Fo] that such examples can not exist). 

The construction of curvature-type invariants whose vanishing is equivalent to the exis- 
tence of a cubic integral is discussed by Kruglikov in |Krl §12]. In particular, he has shown 
that the dimension of the space of cubic integrals is ^ 10, moreover, the metrics admit- 
ting 10- dimensional space of cubic integrals have constant curvature. It was conjectured by 
Kruglikov in jKrj that the next largest value of the dimension of the space of cubic integrals 
is 4 and that in this case g is Darboux superintegrable (the definition see in |K-K-M-W] or 
in [Br-Ma-Ma] §2.2.4]). If the Kruglikov conjecture is true, |Krl Theorem 2] gives curvature- 
type invariants whose vanishing is equivalent to the existence of the 4-dimensional space of 
cubic integrals. 

In the present paper, we answer the question whether the pair (metric g, holomorphic 3- 
codifferential A) is compatible in the sense there exists a cubic integral for g whose Birkhoff- 
Kolokoltsov form is A. On the level on PDE, we reduce the number of unknown functions, 
which of course make the problem easier. This is the reason that our answer (for cubic 
integrals) is more simpler than the answers of [Kr] and [Br-Du-Ea] (for quadratic integrals). 
From other side, the assumption that A is given is natural from the viewpoint of classical 
mechanics. Indeed, all integrals corresponding to different values of h (see §1.41) have the 
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same holomorphic 3-codifferential A. Moreover, on closed surfaces, the space of holomorphic 
3- CO differentials is finite-dimensional. That means that, in the investigation of cubic integrals 
on closed surfaces, one can view our compatibility conditions as the conditions on the metric 
only depending on finitely many parameters. 

2. Complex calculus and prolongation-projection method. 

2.1. Homogeneous polynomial integrals. Let g he a. metric on a surface S and let 
H = be a Hamiltonian given by the kinetic energy and let F{x,p) be its integral which 
is polynomial of degree ^ d in momenta. Then F{x,p) = Yl'j=o^ji^^P) where each com- 
ponent Fj is homogeneous of degree j in momenta. Direct computation shows that each 
Poisson bracket {H,Fj} is homogeneous of degree j -|- 1 in momenta. The uniqueness of 
the decomposition of a polynomial into homogeneous components implies that the Poisson 
bracket {H, F} vanishes if and only if each bracket {H, Fj} vanishes. By this reason, consid- 
ering polynomial integrals of purely kinetic Hamiltonians H = we can restrict ourselves 
to homogeneous polynomial integrals F. 

2.2. Birkhoff-Kolokoltsov codifferential. Every metric g on a surface S admits locally 
so-called isothermic coordinates {x, y) in which the metric has the form g = X{x, y){dx'^ + dy'^). 
If we orient S by means of these coordinates, then the function z := x + \ y will be a complex 
coordinate for the induced complex structure on S, characterised by the following property: 
The multiplication of 2; by i is the rotation of the chart {x,y) by 90°. The existence of 
isothermic coordinates was discovered by B. Riemann. A detailed proof under very weak 
assumption on the metric (e.g., g is merely continuous) can be found in [Beg . 

In forthcoming calculus we use such a complex coordinate z and its conjugate z instead 
of real coordinates {x,y), since most formulas become simpler and more compact. The 
corresponding dual coordinates are p := ^{p^ — '\Py) and p := ^{px + ipy), the function H 
is expressed as H = and the canonical symplectic form in the cotangent bundle T*S 
( "phase space" ) is expressed as Ucan = dp A dz + dp A dz. 

Now assume that F = F{z,z,p,p) is a (local) integral of H cubic in momenta. Then 
F = ^[a-p^ + b-p'^p) where 3? denotes the real part of a complex expression and a = a{z, z), 
b = b{z, z) are some smooth (local) complex-valued functions on 5*. An explicit calculation 
(see |Du-Ma-Toj ) yields: 

{F, H} = ^3f?(/-aj + p^pibXs + 3a\, + \b-, + \a,) + p^p^-{2b\, + \b,)) , (2.1) 

where {■)z and (■)j denote the complex derivatives: 



Since the term ■ must vanish, we obtain immediately the following 

Lemma 2.1. In the situation above, if F = ^{a-p^ + b-p'^p) is an integral, then the coefficient 
a is a holomorphic. 

Rewrite any function F{x^, x^;pi,p2) cubic in momenta as a polynomial F = F'^^^piPjPk in 
Pi. Then, under change of coordinates, its coefficients F'''^^ [x^ ^ x^) transform as components 
of a symmetric (3, 0)-tensor. The same holds for the "a-component" of F: Under complex 
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( -v^ holomorphic) coordinate change z z' it transforms as a' = a ■ (^) . Indeed, the 
momenta p and p are transformed a.s p = p'^, p = p'^ implying 

F = ^ia-p' + b-p^p) = ^(^a-{%Y-p" + b-{^Y %-p'^p'y 

a' 

This means that A := a{z) {-§^)^'^ is a well-defined section of the bundle (T'^^'^-'S')®^ in- 
dependent of the choice of a local holomorphic coordinate. By Lemma 12.11 the section is 
holomorphic. 

Definition 2.1. For any function F{x^,x'^]pi,p2) cubic in momenta the complex tensor 
A := a ■ called the Birkhoff-Kolokoltsov 3-codifFerential associated with F. We shall 

usually denote by a = a{x,y) = a{z,z) its complex-valued coefficient (wrt. some complex 
coordinate z), and denote by A := ^{A) its real part. The latter is a symmetric (3, 0)-tensor. 

Since every non-vanishing holomorphic function f{z) in one variable has isolated zeroes, 
the same is true for any holomorphic 3-codifferential. Notice that if a{z) is holomorphic 
and non-vanishing, then locally near a given point p there exists the root a{z) which 
is again a non-vanishing holomorphic function. It follows that making the change to the 
new complex coordinate z' given by z' := J -ij-j-^ ^ ^^e component a{z) transforms into the 

constant a'{z') = 1. 

Another applications of holomorphicity of A are as follows: 

Corollary 2.1 ( [Kolj ) . i) No Riemannian metric on a closed oriented surface S of genus 
g ^ 2 admits a global non-trivial polynomial integral F . 

a) Let S be the 2-torus T^, g a metric, z a linear holomorphic coordinate, and F = 
3?(a ■ p^ + b ■ p^p) a cubic integral. Then a is constant. 



Proof, i) As is is shown in § 12.11 we may assume that our invariant F is homogeneous 
of degree k. Evidently one can take F in the form F{z, z,p,p) = 9ft(a ■ p*^) + pp ■ B where 
B{z, z,p,p) is homogeneous in momenta of degree k — 2. An explicit computation shows that 
similarly to the cubic case {H, F} = -^^{az ■ + pp ■ C for some function C{z, z,p,p) 
which is homogeneous polynomial of degree /c — 1 in momenta. Thus must vanish which 
means that a{z) is holomorphic as in the cubic case. Similarly to the cubic case, A := a- 
is a well-defined holomorphic section of the bundle (T^^'^^S)'^''. 

The classical complex geometry (see e.g., jGr-Ha] ) says that for any non-zero meromorphic 
section A of the bundle (T^^'^'-'S')'^'^ the difference n{A) —p{A) between the numbers of zeroes 
and poles of A is k{2 — 2g). Since A is holomorphic in our case, there are no poles and 
k{2 — 2g) must be non-negative. Consequently, in the case g ^ 2 A must vanish identically. 
It hollows that \B = F ■ is a global homogeneous integral for H of degree k — 2. Now 
we can apply induction in k. 

ii) A local holomorphic coordinate on a torus is linear if it comes from the standard 
complex coordinate z on the complex plane C under some conformal isomorphism = C/A 
for some 2-lattice A C C. An equivalent condition is that ^ is a globally defined holomorphic 
vector field. Such a field on a torus is non-vanishing, and therefore a is a globally defined 
holomorphic function on a closed Riemann surface. Such a function must be constant by 
Liouville's theorem. □ 
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2.3. Prolongation-projection method. The "naive" idea behind this method is very sim- 
ple. Suppose we have a system S = {Si}, of PDEs which is overdetermined, i.e., the number 
of equations is bigger than the number of unknown functions. In our case, we have two equa- 
tions on one unknown function. We prolong the system, i.e., we differentiate the equations 
with respect to all variables (in our case, x and y) and add the results to the system S. In 
other words, we consider the systems 

S^'^ ■■= {^.} [J{£^^} [J{^y^^} first prolongation, 

S^'^ ■■= {^^} U{|^.} [J&} Ulafe^.} U{^^.} - second prolongation, 

and so on. Obviously, the operation "prolongation" does not change the set of sufficiently 
smooth solutions of the system. Indeed, every (sufficiently smooth) solution of S is also 
a solution of the prolonged systems. From the other side, every solution of the prolonged 
system is evidently a solution of S, because 5 is a part of the prolonged system. 

If the system is overdetermined, then generically the number of new equations grows faster 
than the number of derivatives of the unknown functions. Then, in the generic case, after 
say / prolongations one can resolve the prolonged system iS*^'^ with respect to the highest 
derivatives. In this case, the system S is called the systems of finite type. In our case, the 
system is indeed of finite type, and the number of needed prolongations is / = 1. 

After this point, we can start the prolongation-pro/'ect/on procedure. The key observation 
here is as follows: The successive prolongation S^'"'^^^ can be still resolved with respect to 
the highest derivatives, the number of the equations in iS^'"*"^^ will be generically greater 
than the number of highest derivatives, and with the help of algebraic manipulations we 
can obtain new equations of lower degree. This algebraic manipulation — expressing of 
higher derivatives via lower ones using part of equations and substituting these expression 
in remaining equations — is called projection procedure. We can repeat the prolongation- 
projection many times (clearly, is sufficient to differentiate only newly obtained equations). 

In the generic case, after finite number of prolongation-projections, we obtain a system in 
which the number of algebraically independent equations is greater than the number of partial 
derivatives of unknown functions (we consider the unknown functions as partial derivatives 
of order 0). Such an algebraic system is inconsistent. In this case the initial system S has 
no solution. It may also happen that certain prolongation S^''\ considered as a system of 
algebraic equations on the derivatives of the unknown functions, is algebraically inconsistent, 
even if the number of algebraically independent is less than the number of partial derivatives. 

Thus, the existence of a solution of an overdetermined system S of finite type implies that 
after certain number of prolongations and prolongation-projections we come to the point 
where prolongation-projections do not produce essentially new equations. The latter means 
that the equations we obtain are algebraic corollaries of the equations we already have. 
Moreover, the system of algebraic equations on the partial derivatives of unknown functions 
is consistent. In this case, the system is called involutive, or in involution. One can show that 
(under certain additional assumptions which are fulfilled in our case) involutive systems can 
be solved locally. Moreover, one can find all solutions with the help of solving of certain 
ODEs and of algebraic operations. 

As we noted above, a generic overdetermined system of finite type is inconsistent and 
can not be put to be in involution by prolongation-projections. One obtains a system in 
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involution after certain number of prolongations and prolongation-projections, if the coeffi- 
cients of the initial system satisfies certain partial differential equations (called integrability 
conditions). These PDE on the coefficients are equivalent to two conditions: 

(1) we obtain no new equations after certain prolongation- projection. 

(2) All systems we obtain by prolongations and prolongation-projections, considered as 
algebraic systems on the partial derivatives of the unknown functions, are algebraic 
consistent. 

In our paper, we find explicitly all integrability conditions for the system that corresponds 
that g and A are compatible: they are Q2, G3, ^2 '^3' l^y, Dq from Theorems \T?2\ and 
Propositions EH m 

The prolongation-projection method is a very powerful method for finding solutions of 
overdetermined systems of PDE of finite type. Moreover, in our case the initial system of 
PDE is an inhomogeneous linear system, implying that standard difficulties due to impossi- 
bility of solving explicitly systems of algebraic equations do not appear, so the method can be 
applied algorithmically. Actually, the only difficulty in applying the prolongation-projection 
method for linear overdetermined systems of PDE of finite type is the calculational one. 
In many cases, even modern computer algebra programs are not powerful enough for using 
prolongation-projection. 

We overcome this difficulty with the help of "advanced complex calculus" introduced in 
§2.41 roughly speaking, this calculus is a collection of quite nontrivial tricks that allows us to 
'hide' the covariant derivatives of the objects, and, therefore, makes all formulas very com- 
pact, so we could do prolongation-projections "by hands", without using computer algebra. 
Moreover, all differential equations on the coefficient of the metric and of the Birkhoff- 
Kolokoltsov form are automatically expressed in the invariant form, so one can calculate 
them in an arbitrary coordinate system. 

Notation. Applying the prolongation-projection method we use the following terminology. 
Let £^ be a (differential) equation or expression and S = {Si} a system (ie., a set) of such 
equation. Then S is an algebraic consequence of S if it can be deduced from S by means of 
algebraic manipulations, and a differential consequence if S is an algebraic consequence of S 
and partial derivatives of equations from S. In this paper we work with linear PDEs and 
make linear-algebraic manipulation with equations. In particular, we have linear-algebraic 
consequences of linear systems of equations, and differential consequences can be obtained 
applying linear differential operators. 

Remark. The name "prolongation-projection" is natural in the context of jet bundle ge- 



ometry, we recommend the textbook Kr-Ly-Vi as a source. Equivalent approach is called 



Cartan, or Cartan-Kahler theory, we recommend the textbooks |B-C-G-G-Gt Ilv-La] for more 
details. The system of PDE we consider in our paper is easier than the generic systems 
treated in Kr-Ly-Vi IB-C-G-G-Gt Ilv-La] (because they are linear and of finite type), that's 



why we can achieve our goals without going too far into the geometric theory of PDE. Ac- 
tually, our intention is to make the paper understandable for a standard mathematician 
working in mathematical physics; that's why we avoided the terminology of the geometric 
theory of PDE. 
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2.4. Advanced complex calculus. Our next goal is to develop certain calculus which will 
allow us to simplify computation and express results and formulas in more compact form. 

Let (5", g) be a Riemannian surface with a fixed orientation and z = x + \ y a local complex 
coordinate. Let T*S''^ be the complexified cotangent bundle, its sections are complex 1-forms 
on S. Then sections dz = dx + \dy and dz = dx — \dy form a local frame of T*S'^. Consider 
the decomposition T*S'^ = ^l^^''^^S © ^l^'^'^^S in which the summands are generated by the 
forms dz and respectively dz. The formula for transformations of the forms dz and dz under 
holomorphic coordinate changes shows that the decomposition T*S'^ = il^^'^^S © fi'^'^'^-'S' is 
independent of the choice of a local complex coordinate. A similar decomposition TS"^ = 
7^(1=0)5' y(o,i)5' ^Qj- ^j-^g complexified tangent bundle is obtained using local vector fields 
^ •= ~ '^0_and -§= := + 1^). 

We refer to [Gr-Hal IForj for basic properties of these decompositions. Main properties of 
the bundles T^^'O)^, T^^'^)^, Q^^'^'^S, Q^^'^'^S which we shall use are as follows: 

(1) The bundles n^^^^^S are fi^^'^^S are complex dual line bundles to T^^''^^S and T^^'^^S, 
respectively. 

Using this fact, we define complex line bundles T^P'^^S and ^'^''^IS' with integer p, q setting 

These bundles should be not confused with bundles of tensors of type {p, q) used in differential 
geometry, and with bundles of (p, g)-forms used in complex analysis and geometry. In order 
to emphasise the difference, we use brackets [p,q] instead of parentheses (p, g). 

(1') From the definition and property (1) we obtain immediately 

The bundle T^^'^l = Q'^^'^^ is trivial and its sections are usual complex-valued functions. 

(2) The bundle (T*5)(^'=)©C (respectively {TSY^'^^^C) splits into the sum of line bundles 
each isomorphic to some ^1'^P''^^S (respectively T^''^^S) with p,q ^ and p + q = k. A 
similar splitting holds for a general complex tensor bundle (TS)^^^^ © (2"*5')(®0 ^ 
The components of (T*S)^'^^^ i^C can be indexed by the order of factors il'^'^'^^S 
in the product. 

Notice that the subbundle of symmetric tensors of (J^*S)^'^^^ © C is isomorphic to 
the sum ©p^o^'^''^"^^'^', so that each ^I^P'''~p^S appears exactly once. In the special case 
k = 2 the power (T*S)^'^'^^ © C has as summands one bundle fi'^'^'l, one bundle ^2'°'^!, 
and two bundles Q^^'^^S © Q^^^^^S, Q^^'^'^S © Q^^^^^S each isomorphic to fit^'^l 

Moreover, for positive i ^ min(p, k),j ^ min(g,/), the natural isomorphism 

T^P'^^S © Q^'^'^^S ^ T^p-i'^-^]s © (2.2) 

coincide with the (2+j)-fold index contraction in the space of symmetric complex- valued 
{p + q,k + /)-tensors. 

If 2; = X -|- i ?/ is a local complex coordinate on S, then 

dz ® dz = {dx^ + dy"^) + i {dx ® dy — dy ® dx) = {dx^ + dy^) + \dx A dy 

and dz ® dz = {dx"^ + dy"^) — \dx A dy. Thus the symmetric part of the sum il^^'^^S © 
n(°'^)S'©n(°'^)5'©n(^'°)5' is spanned by dz»dz := \{dz®dz + dz®dz) = dx'^+dy^ and its 
antisymmetric part spanned hj dz Adz = 2'\dx A dy. This gives another decomposition 
of sum n^^''^^S © fi^^'^^S" © (^(O'^^S' © n^^'^'^S into two bundles each isomorphic to l^l^'^l 
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The ant i- symmetric part is the bundle A'^S of complex 2-forms on S and the power 
(T*5)(®2) c splits into the sum of fil^-o]^ ^[0,2]^ 

In what follows we shall use only symmetric tensors. Thus the product dz dz will 
be understood as the symmetric one. This corresponds to the embedding of fi'^'^^ into 
T*S ®T*S (E) C with symmetric image. 

The complex conjugation in the bundles (TS)'^ and {T*S)'^ induces conjugacy bundle 
homomorphisms a : T^^'^'^S T^O'^^^ and a : Q^^^^^S fi^^'^^S with = Id. 

This a extends to conjugacy homomorphisms a : T^''''''^S T^'^'^^S' and a : il'^''''''^S — *• 
fl^i'P^S for any pair of integers p, q. In the case conjugacy a coincides 

with usual complex conjugation. We use the notation a for the complex conjugate of 
any section a. 

If 2; = X + is a local complex coordinate on S as above, then dz^ dz'^ is a local 
trivialisation of rit^'^lS*, and its complex conjugate is dz'^ dz^. 

In the special case p = g = 1 we obtain that dz dz = dx^ + dy^ is real. Further, the 
local expression of the metric in the local coordinates x,y is g = X{dx^ + dy"^). So we 
can consider g as a global section of r^^^'^lS*. 

The metric volume form u = Xdx A dy is the real section of the bundle A'^S ® C of 
complex 2-forms on S. Moreover, under natural isomorphisms fit^'^JS* A^S* ® C the 
metric form g = X^dx"^ + dy"^) is mapped onto the volume form uo = Xdx A dy. The 
difference between r^'^'^lS* and A'^S ® C lies in that how these bundles are embedded in 
T*S(ET*S®C. 

Let V be the metric covariant derivative on 5*. Then it can be applied to any section 
a of any bundle n^'^^S (or T^p^^^S) and the resuh lies in the product Q^p^'^^S » {T*S)^ 
(or respectively T^p^^^S ® {T*S)^) which is the sum Q^p+^^'^^S © Q^p^^+^^S (respectively 
TiP-^^i^S © T^P'i-^^). We denote by V(^'°) and V^^'^^ the component of Va lying in 
Q\p+h<}]S and n^P'^+'^^S, respectively. In "dual" notation, V^^'") and V^"'^) act from 
Tlf'^l^ to TtP-i'^^S and T^p^i-^\ respectively. 

Exphcit calculation show that in coordinates z = x + \y with g = X{dx^ + dy"^) one 
has 

V^^^''\fdzPdz'i) ={§-J -p§-j)dzP+^dz^ 
V^°'^\fdzPdz'') = {-§-. f - q§--f)dzPdzi+^ 

In particular, V'-^'^-'f? = V*-^''^^^' = 0, V*-*^'^-* coincides with ^ for sections of holomor- 
phic bundles T^P'^^S, n^P'^^S, so that such a section a is holomorphic iff V*-^'°^a = 0. 
The operators V*-^'*^-* and V*-^'^-' of complex and anti-complex covariant differentiation 
do not compute. We have 

(i,o)y (0,1) _ ^(o,i)^{i,o))^^^p^-, ^ . ^ . ^ . f^^p^^, (2.3) 

where R is the Gauss curvature of the metric g. The formula is deduced as follows. 
First, one computes the formula directly for the case / = 1. Then using Leibniz rule 
one shows that (V'-^'^-' V*-'^'"'^' — V*-'^'^-' V*-^'^'') is not a differential operator but a bundle 
homomorphism, and hence the formula holds for arbitrary /. 
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(7) To simplify notation, we denote the complex and anti-complex covariant differentiations 
V'-^'^'' and V''°'^^ by {■)-z and respectively. Iterated differentiations are denoted like 
{■);zzzz and so on, the order of symbols z/z coincides with the order of corresponding 
differentiations. In particular, the formula (12.31) now reads 

ifdz^dz'^),,, - UdzMz'i\-zz = ^ -R-g- fdz^dz^; 

^ _ (2.4) 

ifmy,z-z - um);-zz = V ■ ^ ■ ^ ■ 

(8) We shall use the formula 

f;zh;-z- UK = b{f^h}g. (2.5) 

for the Poisson bracket of functions / and h on S with respect to the metric symplectic 
form ujg = \dx A dy = jXdz A dz. Its calculation is as follows: 

Uh,.-U^z = §^Jzd-z=-^^l^\dzd-z=-^g{f^h}g. 

(9) Sections of bundles Vl^^'^^^S, especially holomorphic ones, are called k-differentials. For 
this reason we call sections of bundles T^^'^^S k-codijferentials. 

Let us now recalculate the Poisson bracket {F,H}. Recall that we consider the functions 
F on T*S which are polynomial in momenta. In local coordinates z,z,p,p, such a function 
is a polynomial in p,p whose coefficients are smooth functions on S. With every such 
polynomial Ylij=o fvP'^P' associate a section of the bundle sum 0f j=o -^'^'"'''^ using the 
rule p I— s> -§^,P ^ -§= and extending it in the obvious way on polynomials. Let us denote this 
map by T : F T(F). In particular, T{H) = \g'^. 

Lemma 2.2. T{{F,H}) = T{H) ■ V(T(F)). 

Proof. Let us observe that since both sides satisfy the Leibniz rule in F and commute with 
complex conjugation, it is sufficient to check the formula in the case F = p. The rest 
follows. □ 

In view of the lemma, we can replace the functions on T*S polynomial in momenta by 
their T-images which are finite sums of sections of T'^'-'lS'. In our special case first integrals 
cubic in momenta correspond to sums ^{A + A + B + B) = ^{A + B) where A is a section 
of T^^'^^S and i? is a section of T^^'^^S. The equation 'St{SJ A + VB) = takes values in 
the sum (B'j^^iT^'^'-^'-'^S, and homogeneous components of the equation are: V'-^'^M = 0, 
Y(i.o)^_l_y(o,i)^ = 0, the complex conjugates of these two equations, and 'St{V^^'^^ B) = 0. As 
we have seen above, the equation V'-'^'^M = on a section ofTl^'Ol^ means its holomorphicity. 

2.5. Component B and the principle equation. Let us now consider the last equation 
sfj(y(i:0)^-j = 0. The tensor 5 is a section of T^'^'^^S, so ■ B is section of the bundle 
f2[o,i]5'_ Hence g^ ■ B has the form g"^ ■ B = P(x,y) ■ dz for some complex function /3(x, y) = 
Pi{x, y) + i ■(32{x, y) with real and imaginary components /32- By property ([5]) on pagefT^ 

^2y{l,0)5 = V^m^g^B) = %dzdz. 

Since dzdz is real, the equation 3ft(V^^'°^-B) = is equivalent to 
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This equation can be seen as the closedness of the 1-form — f32dx + (3idy . Hence the equation 
s)fj(^y(i.o)^^ = is equivalent to a local existence of some real function K{x,y) such that 
Pi = ^ and /32 = The latter two conditions can be rewritten as a single complex 

equation /3 = — 2i^, or equivalently 

B = -2\g-^K.,. (2.6) 

The latter equation is identical with the equations (2.7), (2.8) in |Du-Ma-Tb] . 

Substituting (12. 6p in the remaining equation V^^'°M + V*-°'^^-B = 0, we obtain the following 
necessary and sufficient condition: 

Lemma 2.3. Given a metric g and a holomorphic 3-codiff'erential A on a surface S, A is 
compatible with g a if and only if the equation 

S,: K,,, = -y^A,, (2.7) 

has a smooth real-valued solution K, and in this case F = ^{A + B) with B = —2\g~^K.z 
is a cubic integral with the Birkhoff-Kolokoltsov tensor A. 

We call (12. 7p the principle equation. 



3. Proof of Theorem 1.1. 

3.1. Calculation of further equations. We apply the prolongation-projection method to 
the equation (12. 7p considering the function K as an unknown function and g, A (and their 
derivatives) as known parameters. The equations obtained by prolongation-projection pro- 
cedure are called deduced. Instead of considering the real and imaginary parts of equations, 
we shall mostly use complex equations and their complex conjugates. In particular, the 
complex conjugate to (12. 7p is 

K.^zz = ig^A.-, (3.1) 

Differentiating (12. 7p and (13. ip with respect to z and z (first prolongation), we obtain the 
following 4 equations: 

K-zzz = ^g'^A-zz K.^zzz = ^g'^A-zz K.2ZZ = ~^9'^A.zz K-zzz = ^^g'^A-zz (3-2) 

This gives us the expression of all partial derivatives of K of order 3 via lower order deriva- 
tives, which are "hidden" due to the covariant form of the equations and do not appear 
explicitly. In particular, the system {£z,Sz) has finite type. 

Differentiating once more (second prolongation), we obtain 6 equations of order 4: 3 as 
order-2 derivatives of (12. 7p and 3 more as derivatives of its complex conjugates: 

jy- — Lr?A- K -— Ln'^A-- K —— Ln'^A — 

^^\zzzz rftj ■n.-zzz ^^\zzzz 9i/ ^;22:2 ^^\zzzz oi/ ^:2zz . . 

z ' z < z ' f3 3) 

__ — _Ln'^ A K — —1^2/1 _ jy- __i^2/l __ \ - ) 

^^■,zzzz 2" -'^;zzz ^^\zzzz 2" -'^;zzz -'■'■■,zzzz 2" -'^]zzz 

Since there are 5 partial derivatives of K of order 4, at this prolongation-projection step we 
can obtain one new equation. 

Comparing r2['''''']-types of expressions in (13. 3p we see that the desired equation should be 
the difference of equations K-zzzz = ^g'^A-zzz and K.zzzz = -^g'^A-zzz- 

{K-^zzzz ~ K-zzzz) ~ ^g'^{A-^zzz + A-zzz) (3-4) 
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Computing the difference K.^zsz — K.^zzz we use the formula for commutator of covariant 
derivatives, see property ([6]) on page [HI This gives 
K K— = 

^^■,zzzz ^^\zzzz 

{^■^■,zzzz -^;zzzz)l ~l~ {^■^■,zzzz ■,zzzz')2 ~l~ i^-^izzzz -^-jZzzz^S ~l~ i^-^izzzz K -^^zzz^ i 

{{\RgK,,),,)^ - iilRgK,,),,), = ^{K,zR -z - K.-.R,,). 

Here we use fl2.3p which shows that the differences (. . .)2 and (. . .)3 vanish and gives the 
formulas above for the differences (• • .)i and (. . .)4- So using (12. 5p we can write down the 
first deduced equation (new equation in the terminology of 



{K,R} = m{A,,,,)=:D^ (3.5) 

where we denote the right hand side by Dq =: i^^A.^zz) = 2^4. + 2A.2zz- 

Remark. We see that the first deduced equation does not contain second derivatives of K. 
This is not a coincidence but follows from the condition that the higher coefficient in the 



(12. 7p is constant, see |Kr-Ly| for details. 



At this point we can solve the problem of existence of cubic integrals with a given Birkhoff- 
Kolokoltsov tensor A for metrics with constant curvature. 

Proposition 3.1. Let g be a metric metric on a surface S of constant curvature and A a 
holomorphic 3-codifferential. Then A is compatible with g if and only if A satisfy the PDE 
Do = 0. 

Moreover, for a given point P G 5* there exists a unique cubic integral of the form F = 
^{A) + 2H ■ b^pi with a vector field b := b^d^i having prescribed values b{P) = v E TpS and 
rot b{P) = weRatthe point P. 

Proof. Let (a;^,a:^) be any given coordinate system on S. The formula (15. 7p gives the co- 
variant tensor form of the principle equation (12.70 . Denote by Sij the components of this 
tensor form. Let S be the system consisting of equations Sij and their 1st order deriva- 
tives Sij-k- Introduce new variables Kj = Kj(a;^, x^), Kj^ = Kij{x^,x'^), substitute the latter 
into equations Sij instead of the corresponding derivatives K-ij, and denote the obtained 
inhomogeneous linear-algebraic equations by Sij(K). Consider the system S(K) consisting 
of equations £ij(K.), their partial derivatives £ij(K.).k, the equations (Ki)-j = Kjj and the 
consistency equations K12 = K21 and (Kjj).fc = (Kjfc).j. Then the systems S and 5(K) are 
equivalent: the inverse substitution of K.i,K-ij instead of Kj,Kjj transforms S(K) into S 
(turning consistency equations into identities). Moreover, both system are involutive simul- 
taneously. 

Now, since the system S(K) consists of 1st order linear PDEs and is involutive, it admits a 
local solution by the classical Frobenius theorem (see eg. jIv-La] . § 1.9). Moreover, for every 
P E S and every values Kj(P),Kjj(P) satisfying algebraic equations Sij{K.) at the point 
P there exists a local solution (Kj(x), Kjj(x)). For this K-solution there exists a function 
K[x), unique up to adding a constant, which satisfies the relations K-i = Kj, K.y = Kj^. 

To calculate the number of parameters, we use linear-algebraic equations Sij(K) and ex- 
press the functions K12 = K21 and K22 via Ki,K2, and Kn. Substituting these formulas 
in the remaining PDEs, we reduce the number of unknown functions to three. Let S'(K) be 
the obtained system. Notice that some of the substitutions above mean the application of 
prolongation-projection procedure, for example, substitution of the formulas diK^jk = J^ijk 
in equation (Kjj).^ = (Kjfc).j. However, the involutivity of S(K.) implies that every equation 
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in S'(K) is still a linear- algebraic consequence of the system 5(K). In turn, this implies 
the involutivity of the system iS'(K): Indeed, every differential consequence of S'(K) is a 
differential consequence of 5(K), and hence is a linear-algebraic consequence. 

By the construction, the system iS(K) is linear- algebraic equivalent to a system containing 
the equations of the form 9jKj = jFj^, diKjk = J-'ijk where J^ij, Tijk are linear inhomogeneous 
expressions in Kj, Kjj. It follows that the system S'{K) is also linear-algebraic equivalent to 
a system consisting of the equations of the form diKj = JF^,9jKll = JF'^^^ where Tip Tin 
are as above linear inhomogeneous in Kj,Kii, and we may assume that the system 5'(K) 
is of this form. 

Now the generic solution of the system S'{K) is constructed as follows: Given a point 
P and values Ki(P), Ki2(P), Kii(P) at P, we fix local coordinates x := and y := x^, 
integrate the ODE system ((9xKj = Tu, <9xKii = Tm along some interval on the x-axis, and 
then integrate the ODE systems {dyK^i = JF2j,9yKii = ^211 along every interval parallel to 
y-axis. 

It remains to show that the prescribed values as in the assertion of the proposition can 
be used to parametrised a general solution. Indeed, as such parameters we can use the 
values of the 1st order derivatives of K and the Laplacian AgK (instead of K-n at a given 
point P E S. Now recall that by § 12.51 6* and K are related as K-i = uJijV so that h is the 
Hamiltonian vector field on S of the (Hamiltonian) function K (with respect to the form 
ujg). Consequently, h is the skew-gradient of K and the Laplacian ^gK is the rotor rot 6. □ 

3.2. Proof of Theorem 1.1. Let us now turn to the (most interesting) general case when 
R is non-constant. Write the equation (13. 5p as K-^Pz — K-^Rz = ^qDq. Differentiating it 
in z and z we obtain two more equations of order 2 in K: 

K,,,P,, + K,,P,,, - K,,,P,, - K,,P,,, = {Do)z 

K.^-zR-z + K.,R;2, - K,,,P,, - K,,R,,, = ^g-iDo), 

Together with (12.71) and its conjugate we now have 4 equations of second order on K. So we 
expect to obtain at this step one new equation of order 1 on K. It is obtained as follows: 
Add equations (13. 6p with coefficients R-z and R-z'- 

K-zzR;zR;z + K.zR;zzR;z " K.zR;zzR;z + 

K,zR-z-zR;z - K,-z-zR,zR,z - K,-zR,z-zR,z = ^ ' {D o) .zR^Z + yiDo);-zR;Z 

Rearranging it, we obtain 

K,z{R,zR;z):z " K,-z{R.zR -z);z = ^giDo)-zRrz + ^g-{Do)--zR-z + K,-z-zR,zR,z - K,zzR -zR -z. 

The expression R-zR]z equals ^glVPl"^. Let us denote ||VPp by (pi. So the left hand side of 
the latter equation is ^g^{K, (pi}. On the right hand side we have ^g[{Do).zR-z + {Dq).zR-z) 
which equals ^g^{'VDo,'VR). In the last two terms we use (12. 7p . So finally we obtain the 
equation 

{K, (^1} = (VZ^o, VP) - m{A.,z ■ {R,zf) =: Z^i, (3.7) 
in which we denoted the right hand side by Di. 

From this points we can produce many equations, each of the form {K, (p} = D where 
ip and D are real functions which are certain differential expressions involving the metric g 
and the tensor A-z- In fact, we have two such "dummy" procedures: 
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The first one is the formal repetition of the deduction fl3.5p ^^ fl3.7p and uses the fact that 
the exphcit expressions for R and Dq were not involved. So from any equation {K^ ip} = D we 
can obtain a new equation {/T, (yj'} = D'wiih.Lp' := i|Vv3pandD' := {V D,V'p)—4Si{A.z'f%)■ 
T\le second procedure is based on the Jacobi identity for Poisson brackets: starting with 
any two equations of the form {K, ^p'} = D' and {K, p"} = D", we obtain 

{K, W, p"}} = {{K, p'}, p"} - {{K, p"}, p'} = {D', p"} - {D", p'} (3.8) 

and we can set p'" := {cp', cp"} and D'" := {D', p"} - {D", p'}. 

Let us notice that both procedures are application of prolongation-projection Method: In 
the first case we differentiate the equation {K, p} = D in z and z, add the equation (12. 7p and 
its complex conjugate, and make linear-algebraic manipulations on the system of 4 equations 
excluding second order derivatives of K. In the second case we differentiate both equations 
{K, p'} = D' and {K, p"} = D" and then make the same linear-algebraic manipulations on 
the system of 4 equations. 

To unify the notation, we set 

Pq:=R, := i|Vv2oP, p2 ■= {^q.P^i}, p^:=\\Vpi\^, 

Do := 4 3fJ(A,,,), D, := (VZ^o, V<^o) - 4 3?(A, ■ ((<^o);.)'), (3.9) 

:= {Do,pi} - {D,,po}, D, := (VZ^i, Vy^i) - 4 3? (A, ■ {{pi),,f ). 

Thus the equation {K,p2} = D2 is obtained from {K^po} = Dq and {K,pi} = Di using 
Jacobi identity, and the equation {K,p3} = D3 from {K,pi} = Di using (13.51) and its 
complex conjugate. 

It appears that the system consisting of 4 equations Sp. {K,pi} = Di, i = 0, . . . , 3 is 
involutive and differentially equivalent to the original equation S^- We state a more general 
property which will be used also in the proof of Theorem 11.21 

Lemma 3.1. (a) Let the coefficients of the equations £' = {K,p'} — D' , S" = {K,p"} — D" , 
and £"' = {K, p'"} - D'" are related as p'" := {p', p"} and D" := {D', p"} - {D", p'}. Then 
the equation E'" is a linear-algebraic consequence of the equations £'.^,£'.y, £'.'^,£'.'y. 

(b) Let the coefficients of the equations £' = {K,p} — D and £* = {K,p*} — D* are 
related as p* := l\Vp'\^ and D* := {VD',Vp') - A^{A.,{p'.J^) . Then the equation £* zs a 
linear-algebraic consequence of the equations £!^,£!y, and the equations 

(I) £z and its conjugate £z- 

(c) Under hypotheses of parts (a) and (b) assume additionally that {p', p"} = p'" is non- 
vanishing at a generic point. Let also the equation £'^ = {Kjp''} — D'^ he defined by p'^ : = 
||V(/?"p and D'^ := {V D'\V p") — A'St{A.^z{'^"^Y) . Then the equations (I) are linear-algebraic 
consequences of the following set of equations: 

(II) equations £"',£*, £'^ ; 

(III) the first order derivatives £[^, £[y, £",^, £'.'y. 

Recall that "linear-algebraic consequence" means that new equations appear as linear 
combinations of old ones with coefficients which are rational functions of coefficients of the 
old equations. 
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Proof. Clearly, parts (a) and (b) are simply a restating of the "dummy" procedures intro- 
duced above. 

(c) Rewrite the equation S' and S" in the form Kzip'.^ — K^^'.^ = '2 '\gD' and K^ip'-'^ — Kzip'.'z = 
2\gD". Then the determinant of the linear system is y^'.^y^'-'z ~ v"zf'-z — ^2i(7 ■ {(p',if"}, 

f D' .\ 

and the solution is given by the formulas Ki = J^i with jFj := ■ det( ^„ _Ji 1. 

Differentiating the formulas, we conclude that there exists a unique solutions of the equations 
(111), considered as a system of linear equations on K.^^i^i^ and this solution is given by 
K-^x^xi = {^i);x3- In particular, [J^^.^^ = {J^2);x^- Notice that each (jFj).^j has the form 
^^„p for some PDOs J^ij which are polynomial in (yj', (yj", D' , D" and their derivatives. This 
shows the fact that the determinant of the matrix of leading coefficients of the system (111) 
equals {</?', Lp"Y (up to sign). Indeed, the solutions i^jx^xJ = {^i);x3 can be obtained by means 
of linear-algebraic manipulations with equations (111). 

The key point in the proof is that we can obtain the equations (11) substituting the 
expressions K.^ixi = {^i);xj iu the equations (1). This is also linear-algebraic manipulations 
with equations. The assertion (c) claims that this operation is invertible, which is true. □ 

Define the expressions Qi using the following Jacobi-like expressions: 

Go = {<^i, <^2} ■ ^3 + W2, V^s} ■ Di + {V93, ipi} ■ D2, 

Gi = {Vo, V2} ■ D3 + {v?2, V^s} ■ ^0 + {^3, Vo} ■ D2, 
Q2 = {v^o, V^i} • D3 + {(pi, LPs} ■ Do + {v?3, V^o} ■ ^1, 
= {'Po, Vi} ■ D2 + {(pi, (P2} ■ Do + {ip2, fpo} ■ Di. 

Up to normalisation, the expressions Qi are the (3 x 3)-minors of the extended matrix of the 
coefficients of the equations £0, . . . with z-th row excluded. In particular, 

({.Vo);z {.Vo);z Dq\ 

Q2 ■= 4^^* ('^i)'^ ('^i);^ ^1 ^3 4"^^* 

\(<^3);^ {^3);z D3 J 

Every expression Qi is a PDE on the metric g and the 3-codifferential A. 

Corollary 3.1. Assume that ip2 is n on-vanishing. Then the system S := {So,£i) is in- 
volutive if and only if the integrahility condition Q3 = is fulfilled. A solution K of this 
system satisfies the equation K-zz = —^g'^A-z if one? only if the integrahility condition Q2 = 
is fulfilled. 

At this point we can give the proof of Theorem 11.11 

Proposition 3.2. Let g = A(x, y){dx^+dy^) he a metric and A a holomorphic 3-codifferential 
such that the differential expression ip2 = {ipojipi} is non-vanishing. 

Then g admits a cuhic integral of the form F = ^{A + hp^p) if and only if A satisfies the 
covariant PDEs ^3,^2- Moreover, in this case the complex valued function b = bi + \b2 is 
given hy the formulas bi = A~^/C25^2 = — A~^/Ci with 

{Lpo^fl} \Dl -{Vl);xO 



/ifo);z 


i^o);z 


Do\ 


i¥>l);z 




dA 


\M;z 


M;z 


D2I 
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Proof. By Lemma 12.31 from § 12.51 the existence of a cubic integral F of the form above is 
equivalent to the existence of a solution K of the principle equation K.^g = —^g'^A-z, and in 
this case F = ^{A + B) with B := hp^p = —2\g~^K-z is the desired cubic integral. 

As we have shown above, equations Sq,Si are differential consequences of the real and 
imaginary parts of Sz- So by Corollary 13.11 equations Q2 = 0, ^3 = are necessary and 
sufficient conditions for solvability of Sz- The formulas fl3.10l) give simply the solution of a 
linear system (^o^^i)- D 

Let us observe that under condition of non- vanishing of all brackets {(pi,(pj} any two of 
the equations Qo, . . . ,^3 are linear-algebraic consequence of remaining two. 

4. Metrics with {R, |V-R|^}g = and proof of Theorem 1.2. 

In this section we consider the degenerate case when the function ip2 introduced in (11.61) 
vanishes identically. This will be a general assumption in this section unless the opposite is 
stated explicitly. In this case we additionally suppose that the Gauss curvature i? = (/jq is 
non-degenerate, ie., the gradient of R is non- vanishing. 

4.1. Proof of Theorem 1.2. Recall that the condition = means that the Gauss 
curvature R and the square of its gradient | V-Rp are functionally dependent, i|V-Rp = f{R) 
for some function /(r) of one variable. Notice that the function (p^ = ||V(/?ip is also 
functionally dependent, since 

As in the non-degenerate case, we want to apply the prolongation-projection method. Our 
starting system of equations is again So := {£z,Sz}- Let us analyse the "forking point" in 
the procedure. Since we assume non-vanishing of V-R, there will be no divergence in results 
until Corollary 13. II At this point the degenerate case (f2 = differs from the non-degenerate 
one in two aspects. 

The first one is that we must replace Si by another equation. We shall find such an equation 
in a moment. The second aspect is that in the case (^2 = the result of the projection 
procedure is slightly different as before. Namely, now the system {So,Si) is degenerate at 
every point since the left hand sides of the equations, {K, ipo} and {K, (pi}, are Lie derivatives 
of K along proportional vector fields. Since the equation Si still must be satisfied, we obtain 
a new differential condition on g and A: the right hand sides of Sq, Si must be proportional 
with the same coefficient as the left ones. Clearly, this condition is simply the application of 
the projection procedure to the system {So,Si). 

The resulting equation reads: 

Vz := ((po);.^i - iVl);zDo = 0. (4.1) 

Notice that even if this is formally a complex-valued equation and so two real ones, the real 
and imaginary parts are equivalent under the condition ip2 ^ 0. This means that the system 
of conditions {ip2,T>z} contains only two independent conditions. 
In a more explicit form the equation Vz reads: 

Vz = R.,z ■ (VDo, VR) -R.,z-4- ^{A,z ■ {R;zf) - Do ■ (|| Vi^H^,. 
Another possible form for this condition is the differential 1-form 

det ( ^° ) = Vzdz + V-zdz = V^dx + Vydy, 
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where Dz,T^xj^y are defined in tlie obvious way. One can use any of tlie equation Vr^^Vy 
instead of provided the corresponding derivative ^f- is non-vanishing. 

Now we seek for the substitute for the equation Si. Since our previous step was the 
projection, the next one is the prolongation of the equations. Let S be the system of 
equations obtained so far. Those are Sz,£z and their derivatives up to order 2, the equation 
So and its derivative {So)-x, and the equations ip2 and Vz. 

As the next step we are going to add to S the second order derivatives (£^o);a:*xJ - Here we 
make the following easy observation. It follows from Lemma [XT] and the above consideration 
that the equations {So)-y is a linear-algebraic combination of the equations from S of order 
^ 2 in K and the condition T>z. Similarly, {Sq)-x can be obtained as a linear algebraic 
combination of the condition Vz and the equations from S of order ^2 in K, in which {Sq).^ 
is replaced by {So)-y Thus replacing {So)-x by {So)-y we obtain a system equivalent to S 
provided Vz = 0. 

Consequently, under condition Vz = the derivative {So)-yy is a linear algebraic com- 
bination of the same equations and its derivatives in y. In particular, {S^j-yy is a linear 
algebraic combination of equations from S and the equation {8o)-^xy Interchanging x and 
y we conclude that (£^o);xx is also a linear algebraic combination of equations from S and 
the equation {Eoj-xy In fact, each of the equations {£q)-x^x3 is linear-algebraic equivalent 
to each other modulo the system S and the condition Vz- Thus adding to S an arbitrary 
single derivative (i^o);x*xi instead all three we obtain an equivalent new system. As such a 
derivative we choose {£o)-xy = (^o);zz- 

Now we make explicit calculation. The derivation (£^o);zz gives 

K -jZzzR^z ~l~ K -jZzR^zz ~l~ K -fZzR-^zz ~l~ ^ -jzR-jZZZ /A r\\ 

-K,zzzRz - K,zzR;zz - K,zzR,zz ' K,zR-,zzz = 2\ Q ■ {Do),zz. ^ ' ' 

Simplifying the obtained equation we apply the following relations: the rule (jH]), the ho- 
molorphicity equations A-z = 0, A.^ = 0, and the substitutions fl2.7p . (13.21) . Let us notice 
that the latter case we exclude higher order derivatives of K and so apply the projection 
procedure. Calculating, we obtain K-zzz = ^g'^A-zz and 

K-zzz = K-zzz + {K-zzz ~ K-zzz) = ~^g'^A-zz — ^gRK-z- 
Besides we also use the relations 

R;z-z = R;zz = \9^R aud R-^zzz = R.zzz + \gRR;z = \g^R;z + \gRR;z 

where A = is the metric Laplace operator. Substitution of these relations yields 
y^A--zzR--z + y^A--zR--z-z + \K.zgAR-^-z + 
{^g^A-zz + lgRK.z)R;z + ^^A-zR^zz - K-zHgAR-z + \gRR;z) = {g^ ■ ADq. 
Rearranging the terms and dividing by ^g, we obtain 

K-z/^R-rz - K--zAR.z = 2\gADo -2\g{A.zzR;z + A.zR;zz + A.-zzR,, + 
which finally yields 

{K,AR} = ADo-2^{A,zR,z);z). 
We denote this equation by £1 and set 

:= AR Dl := ADo - 2^{{A,zR;z);z) ■ (4.3) 
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The equation £^ has the already famihar form 

{K,^l} = Dl, (4.4) 
and so we can apply Lemma 13.11 This gives us two more equations 

s; := {K, cp;} - d; s; ■.= {k, cp;} - d; (4.5) 

in which 

:= {^0, ^1} D; := {Do, ^1} - {D^ ^0} (4.6) 

:= i|V^tP Dl := (VZ^I, -4 3fJ(A, ■ ((^D;.)')- (4-7) 

We also define the differential expressions 

^*:= 1 .detff?° "hl^^'Y (4.8) 

Proposition 4.1. Let g = X{x,y){dx'^+dy'^) be a metric and A a holomorphicS-codiff'erential. 
Assume that g satisfies the differential condition (^2 = and that pf^ = {R, AgR}g is non- 
vanishing. 

Then g admits a cubic integral of the form F = ^{A ■ + b ■ p^p) with the given tensor A 
if and only if X and A satisfy the covariant PDEs ^3,^2 , ^'^^ Moreover, the component 
b = bi+ \b2 is given by the formulas bi = A~^/C2 5 ^2 = — A^^/CJ. 



Proof. As above, the existence of such a cubic integral F is equivalent to solvability of the 

ig^A.,. Since V,g*,g* 



equation S^: K.^z = —-kg^A.^. Since V^Q^^Ql are differential consequences oi £z, they are 



necessary conditions. 

Vice versa, the system {£^.,£1) is solvable if and only if the integrability condition ^3 = 
is fulfilled, and then the solution is given by the formulas (14. 8p . Further, by Lemma l3.ll 
(c), the equation £z is a linear algebraic consequence of the conditions of the 1st order 
derivatives (£o);i5 (^r);i and the equations £i,£2,£l. As we have shown, in the presence 
of £o,£i those three equations are equivalent to V,Q^,Q2, respectively. Thus under the 
conditions V = 0,^2 = ^^^3 = the solution K of the system {£0, £^) solves also the 
equation K.^z = —■^g'^A-z- The proposition follows. □ 

4.2. Metrics admitting a Killing vector. In previous paragraphs we considered the cases 
when one of the differential expressions = {R, || V-R|^}g or ip2 = {R, AgR}g is non-zero. 
Here we treat the problem of detecting of cubic integrals in the case when both ip2 and ip2 
vanish. Recall that by Bonnet-Darboux-Eisenhart theorem (see § 11. 5p this means that the 
metric g admits a Killing vector field L^^, and then L = Upi is a non-trivial linear integral. 

We maintain the notation introduced above. Define the expressions from 
the relations 

In particular, 

V: := {ipo);.Dl - {vl),.Do = R;. ■ {ADo - 2^{{A,zR;z);z)) - AR,^ ■ Do (4.9) 
and similarly for V*,V*,V^. 
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Proposition 4.2. Let g be a Riemannian metric on a surface S with the curvature R and 
A a holomorphic 3-codifferential. Assume that R-x is non-vanishing and that the expressions 
(p2 = {R, ||V-R|^}g and = {R, AgR}g vanish. Then A is compatible with g if and only if 
it satisfies the equations V^. = and V* = 0. 

Proof. Applying the prolongation-projection method, we repeat the proof of Proposition 14.11 
until computation of the equation S^. At this point the hypothesis 2 7^ of the part (c) 
of Lemma 13.11 is not satisfied. The situation here is similar to that in Section H] where the 
hypothesis ip2 was not fulfilled. 

At this step prolongation-projection method produces the equation P* = 0, which is an 
analogue of the equation Vx = in the case = 0. Indeed, the condition (^93 = means 
the degeneration of the matrix of coefficients of the equations So, Si whereas the pair of 
equations {ip2jV*) means the degeneration of the extended matrix of coefficients of So,Sl. 
In particular, under the hypotheses R-x 7^ and = the equations Si = and P* = 
are equivalent. 

Since Vx,T>* are differential consequences of the equations Sz,Sz, their vanishing is a 
necessary condition for solvability of £^2- 

Vice versa, assume that Vx = V* = 0. Let 5* be the set of the following PDEs: Sz,Sz, 
their derivatives up to order 2, the equation Sq and its derivative (^o);x- Then 5* is involutive 
and hence solvable, and every solution K produces a cubic integral. 

The whole set of solutions of 5* can be constructed in the same way as it was done in 
the Proposition 13.11 We fix the initial value K.x{P) at some point P G 5, solve the ODE 
■§-K.^x = J(^fx with the initial value K-x{P) along x-axis, then the ODE ^K,x = IC* with 
the initial value K.x{x, 0) at y = 0, and finally express K.y using K.y = /C*. □ 

5. Invariant expressions. 

The complex calculus introduced and applied in previous sections relies on the choice of 
isothermic coordinates. In this paragraph we get rid of it: Lemma 15.11 gives an answer on 
the following two questions (for a fixed given metric): 

(1) Given an integral F cubic in momenta, how to calculate the tensor A (the real part of 
the Birkhoff-Kolokoltsov form)? 

(2) Given a symmetric (3, 0)— tensor A, how to understand whether it can be the real part 
of a holomorphic 2-codifferential. 

If the coordinates we are working in are isothermal, it is easy to answer both questions 
using the definition. But if the coordinates are generic, the questions are not that trivial. 

We would like to remark that the first question will be especially interesting in view of 
program for searching new metrics admitting cubic we integrals suggest in the conclusion 
Section [71 

We use the following notation, x = and y = x"^ are local coordinates, g = gijdx^dx^ is 
the metric tensor, A := ^ydet{gij) is the volume density, so that Ug = Xdx A dy is the volume 
form, and Jf is the operator of the complex structure, ie., the operator of rotation by 90° in 
the tangent bundle (it is easy to construct it for an arbitrary metric). 
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Lemma 5.1. i) Let F = (F*-'^') be a symmetric 3-vector (ie., {3,0)-tensor). Then it can 
be uniquely decomposed into the sum F = A + B where A = (A'-''^) and 13 = (i?*-'^) are 
also symmetric {3,0)-tensors, and A = ^{A), B = 3?(i?) for some sections A of the bundle 
7^[3,o]^ anc? B of the bundle T^'^'^^S. Moreover, for arbitrary 1-forms q;i,q;2,«3 one has 

A{ai,a2, as) = -(F(ai, 02, as) - F{Jai, Ja2, as) - F(Jai, a2, Jas) - F{ai, Ja2, Jas)) (5.1) 

-B(ai,a2,as) = ^(3F(ai, 02, as) + F(Jai, Ja2, as) + F( Jai, a2, Jas) + F(ai, Ja2, Jas)) (5.2) 
or in the index form 

Q^ik ^ 1 (3^*^fe + F^'^'^'JiJ], + F''^'''ji,4 + F'^''''jj,4,) (5.4) 

a) Let A be a section of the bundle T^^'^^S and A := 'Si{A) its real part. Then the imaginary 
part Q{A) is given by 

C^iA))'^^ = - {A^^^Ji + A^^^Ji + i'-'Vf ). (5.5) 
3 

Further, A is holomorphic if and only if the tensor A = [A^^'') satisfies the equation 
Vv4( J-, J-, J-, J-) = — Vv4(-, ■,■,■) or in the index form 

A^^"^'^ = -A''^''''^'4jp^,jl (5.6) 
where (ijkl) means the symmetrisation in the induces. 

■ 2 

Hi) The principle equation K.^z = —^A.^ can be written as 

IK,m{^^6] + J\fj) + gki^u:j)lA^''^■,m = (5.7) 
where (^ij) means the symmetrisation in indices. 

Proof, i) By property ([2]) from § I2.4[ the bundle of complex-valued symmetric 3-vectors is 
naturally isomorphic to the sum T^^'^^ S (BT^'^''^'^ S ®T^^''^'^ S (BT^'^'^^ S . The complex conjugation 
interchanges two inner and two outer summands. Therefore every real symmetric 3-vectors 
has the form F = 'Si{A + B) for uniquely defined sections A of T^^'^^S and B of T'^'^'S*. 

The formulas for A = ^{A) and B = ^{B) are subject of linear algebra, therefore it is 
sufficient to check them for the Euclidean case (M?, gEuc) and constant tensors A = ^ ■ ^ ■ ^ 
and 5 = f ■ f ■ i. In this case A = 5J(A) = ((|:)3 - 3(|:)^ J), B = ^{B) = I 
(^)^), and J(^) = J{§^) = and an explicit verification follows. 

a) Consider the covariant derivative VA = ^(VA). By ([5]) from § I2.4[ in any complex 
coordinate z we obtain VA = A.^dz + A.^dz. Again in the coordinate z, multiplication 
with g-'^ = A"^^^ gives g'^VA = A^A"^^ + A.z>^-'^^. Thus g~^VA is the sum of two 
components of the type T'^'^'S* and T^^'^^S, in which the T'^'^'lS'-component is "responsible" for 
Cauchy-Riemann term A-z- These two components are easily distinguished by the operator J: 
C{J-, J-,J-, J-) = +C(-, ■', -, ■) for every section C of the bundle T^^'^^S, and C( J-, J-, J-, J-) = 
—C(-, ■,■,■) for every section C of T'^'^'S". Thus holomorphicity of A is equivalent to the 
equation g-'^VA{J-, J-, J-, J-) = -g-'^VA{-, -, -, ■). 

Now we make use two easy observations. The first is that if C is a sum of two sections 
of the bundles T^^'^^^S and T^^'^^S respectively, then the T'^-'^JS'-component of C vanishes if 
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and only if C := 3?C satisfies the equation C{J-, J-, J-, J-) = —C{-, ■, ■, ■)• This fact can be 
deduced from the action of the complex conjugation on bundles T^''^^S, see ([3]) from § 12.41 
Another observation is that in the index form the symmetric tensor ^{g'^V A) is given by 
j^m;i) ^ A^'^'',rng^^"'- This implies the formula CT . 

Now consider the tensor A.^dz. By the construction, this is the section of the bundle 
2^[3,o]^ (g) (^[1^0] 5". xhe latter bundle is isomorphic to T'^'^^^S*, and the isomorphism is induced 
by the isomorphism T^^'^lS* ® r2'^''^^S' = C. The latter isomorphism is simply the operation 
of contraction of indices. This operation commutes with complex conjugation, and therefore 
with the operator 3? of taking real part. It follows that the image of 'Si{A.zdz) in Tt^'^'S" is 
given by the symmetric (2, 0)-tensor divA with components [A^^^'-k). It should be noticed 
that besides the usual symmetry A^^^-k = A^^^.^ this tensor has another symmetric property, 
namely divA(J-, J-) = — divA(-, ■). Thus divA has two independent components, as it should 
be since A-z is complex-valued. 

The last formula expressing '^{A) can be deduced in a similar way as above and is left to 
the reader as an easy exercise (and is not important for us). 

iii) The real part of the principle equation (12.71) is 

K.,,2 + K,z. = ^{A,2-A,,). (5.8) 

Since the principle equation is the f2['''^]-component of (15. 8p . both equations are equivalent. 
Since T{Jv,Jw) = T(v,w) for every section T of the bundle fit^'-'^lS' and T{Jv,Jw) = 
—T{v, w) for sections of r^'^'^lS* © fit^'^'S*, the tensor K-zz + K-^z has components 

Rewriting the right hand side of (15.81) we start with the equality ^{A-z — A-z) = ^(A-z)- 
Since A is holomorphic and A-z = 0, the tensor A-z equals A^^^-^k- Using the fact that the 
isomorphism (12.21) is simply the contraction of indices, we obtain the formula QikQ ji^{A^^^ -^rn) ■ 
Further, since A-^z has type Tl^'Ol^, i-A.^ = J-A.^ and hence 3(A.^) = -3ft(J-A.^). Using the 
equality gijJ^k = ^ik we obtain the desired form (15. 7p of the equation. □ 

6. PSEUDO-RlEMANNIAN CASE. 

In this section we show that the results and formulas obtained in the case of Riemannian 
metrics remain valid in the pseudo- Riemannian case after an appropriate modifications which 
we describe. 

6.1. Null-coordinates. These are the counterpart of isothermic coordinates in the pseudo- 
Riemannian case. 

Definition 6.1. Let g he a. pseudo- Riemannian metric on a surface. A vector (field) v 
satisfying g{v,v) = is called a null-vector (field). Null-coordinate^ are such coordinates in 
which g has the has anti-diagonal form g = X{x, y)dxdy. 

Lemma 6.1 (Folklore). Let g be a pseudo-Riemannian metric on a surface S. Then at each 
point on S there exist local coordinates x = x^,y = in which g has the anti- diagonal form 
g = X{x,y)dxdy 

^From the mathematical point of view, isotropic coordinate system would be probably a better notion. 
However, this terminology is already established in the physical literature, see e.g. |Cu-Laj . and moreover, 
the notion isotropic coordinates itself has another meaning, see e.g. [Cro| . 
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Proof. At each point p on S the metric g has exactly two isotropic directions, ie., pointwisely 
linearly independent vectors fields v,w with g{v,v) = g{w,w) = 0. Integrating them we 
obtain a 2-web of curves on S. The local parameters functions for this 2-web are the desired 
coordinates x = x^,y = x"^. □ 

Let us locally fix the order and the orientations of null-directions. Then an null-coordinate 
system is defined uniquely up to oriented reparametrisations x'^ = 

6.2. Calculus in null-coordinates. As in the Riemannian case, formulas and calculus 
become much simpler if we use null-coordinates. 

In analogy with the complex calculus, define real bundles T^^'^^S, T^^'^^S, Q^^'^^S, 
generated respectively by vectors and the forms dx^, dx"^. The definition of these 

bundles depends only on the local combinatorial data, whereas the decompositions TS = 
Tl^fi]g ^ rp[o,i]g T*S = © n^^'^^S are defined solely by the metric g. Further, for 

p,q ^ set 

Then 

T^P'^^S ® T^p''i'^S = T^P+p''i+i'^S (6.1) 

for non-negative p,p',q,q'. We extend the definition of the bundles T^P''^^S, Q,^P''^^S for all 
integers p, q by means of the formula fl6.ll) and the relation 

qMS = T^~P'-i]s. (6.2) 

Every section of the bundle T^p^^^S or Vi\P''^^S is given by its single component/coefficient 
which is a real function. 

As in the case of complex calculus, we imbed the bundles T^P'^'^S^ VI^-P'^'^S with p, g ^ in the 
bundles of symmetric {k, 0)- and (0, k) tensors. The imbedding is done using symmetrisation 
operator. This gives us natural decompositions Sym'^TS = ®^^g^,^T^P''^^S and Sym^T*5 = 

The metric g is a section of the bundle ri'^'^lS'. We use multiplication with powers of g to 
define natural isomorphisms T^^'^^S ^ T^p+k,Q+k]s ^nd ^[^'"15 ^ Q^p+'^'I^''^ S . 

Further, define the operators V[^'°] : Q^P'^^S ^: Q^p+^'I^S and Vl"'^] : Q^^'^^S ^: Q^P'1+^^ 
as the corresponding homogeneous components of the operator V. Thus for any section 
F = f(x, y)dxPdy'^ of the bundle n'^'^'S" with p, g ^ we obtain 

\/i^fi^F = f,^dxP+^dy\ V'^-'If = f,ydxMy''+\ 

and ^7^^'"^'^ + V'^'^^-F = Sym(VF) where Sym denotes the symmetrisation operator. In the 
case of a section A = a(x, y){dxY{dyY of the bundle T'^'^lS* with p, g ^ we obtain 

and V^^'^^A -|- V'^'^'A = Tr(Vv4) = div(74) where Tr denotes the contraction of indices. The 
operators V'^'*^' and V^^'^^ do not commute, and the commutator formula is similar to the 
complex case (12. 3p : 

2 (6.3) 
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Let F{x,y\PxiPy) be a function on T*S cubic in momenta. Then it has a form F = 
AipI — A2Py + {BiPx + B2Py)pxPy for uniquely defined sections Ai of T'^'^^S', A2 of T^^'^^S', 
Bi of T'^'^'S*, and B2 of T^'^lS". The Hamiltonian function corresponding to the metric has 
the form H = 2\(x\) ' ^"^^ equation of integrabihty of F reads 

{F,H] = ^[pi\- {Ar),y - p4 A . (A2);.+ 

plPy ■ {Bl ■ \,y + 3Al • A;^ + A • + A • (^l);x) + 

(6.4) 

PxP^ • {B2 ■ X-x - 3^2 • A;j, + A • (52);x - A • {A2)-y) + 

pIpI ■ (2^1 • A., + A • + 2B2 -X-y + X- iB2)-y)) 

Lemma 6.2. Let g = Xdxdy he a pseudo-Riemannian metric on a surface S in null- 
coordinates. 

(a) A function a{x,y) is independent of the variable y if and only if for some/every k ^ 
the tensor A := ad^ satisfies the equation A-y = 0. 

(b) Let F be an integral of g cubic in momenta. Then its components Ai and A2 depend 
only on one corresponding variable. 

Moreover, if Ai, A2 are non-zero at a given point on S , then there exists a null- coordinate 
system {x,y), unique up to exchange of the coordinates, in which components of Ai are 
identically 1, ie., Ai = and A2 = dy. 

Here A-y means the covariant derivative V^^'^'A. 

Remark 6.1. Similar statement holds for polynomial integrals of other degree, see |Bo-Ma-PuH 
§3.3] for the case of quadratic integrals. 

Proof. If {x', y') is any null-coordinate system with Ai = ai{x')dl, and A2 = a2{y')dy,, then 
the desired coordinates are defined by x = f „ f^' and v = f , □ 

Definition 6.2. For any function F on the cotangent bundle T*S cubic in momenta we call 
the components Ai and A2 the Birkhoff-Kolokoltsov 3-codifferentials associated with F and 
pseudo-Riemannian metric g. Ai and/or A2 is said to be quasi-holomorphi^ii Ai = ai{x^)d^i 
with aj(x*) depending only on one variable. 

The null-coordinate system (x, y) in which Ai = and A2 = dy is called adapted to 
Ai,A2 or to F. 

Proposition 6.1. Let g = Xdxdy be a pseudo-Riemannian metric on a surface S in null- 
coordinates and F a function cubic in momenta. Then F is an integral for g if and only 

(A) its components Ai,A2 are quasi- holomorphic; 

(B) there exists a function K{x,y) such that the components -81,-82 are related as Bi = 
-X~^K.y, B2 = X~^K,x; 

(K) the function K satisfies the equations 

K-xx = g"^ {A2)-,y K.yy = g^ {Ai).x- 

■^In view of pseztdo- Rieniannian metrics, pseurfo-holomorphic would be probably a better notion. However, 
this terminology is already reserved, see [Gro] . 
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Proof. Five homogeneous components of bi-degree T^^'^^S {i + j = 4) of the Poisson bracket 
(16. 4p are five equations on Ai, Bi and A. The function F is a cubic integral for g = \dxdy iff 
all five equations are satisfied. The first two are equations (A) . The last one can be written 
as d{\'^Bidx + \^B2dy) = and hence is equivalent to the property (B). Finally, substituting 
expressions (B) into the remaining two equations (second and third lines in (16.41) ) we obtain 
equations (K). □ 

6.3. Calculation of equations in pseudo-Riemannian case. Here we use the trick and 
translate all the results of §2.31 and Section H] almost literally in the new situation. Doing so, 
we must follow the next rules. 

• the complex coordinate z is replaced by the null-coordinate x, the coordinate z by y; 

• the holomorphic 3-codifferential A is replaced by the quasi-holomorphic 3-codifferential 
Ai, A by ^2, 

• the imaginary unit i is dropped (replaced by 1), the operator J*j of rotation by 90° is 
replaced by the identity Id operator; 

• the complex conjugation a i— > a is replaced by the null-involution I : {x,y) ^ iy,x); the 
operator Q'(A) = ^(a — a) transforms into the operator : |(a — J(a)), this is the 
operator of the projection on the eigenspace E/(— 1) of the operator / corresponding to 
the eigenvalue —1; similarly, the operator 3fJ(y4) = |(a + a) transforms into the operator 
P+ : a I— s> |(a + /(a)) corresponding to the eigenvalue +1; 

• the Poisson bracket {f,g} = ^{f-zQ-z - f--z9-,z) = g^^'^iUg-,?) is now given by the new 
formula {f,g} := g'^P-{f;x9;y) = j^{f;xg;y - f;y9;x); 

• if it is possible, every equation must be written in complex form and then its conjugate 
must included to the system; exceptions are the equations which are real "in their 
nature", such as Sq, . . . ,Ss;Si, . . . ,S^,Qo, . . . then after applying the above rules 
every pair of complex conjugate equations transforms into a pair of equations S,I{S) 
where / is the null-involution above 

For example, the equation K.^z = A-z is complemented by its conjugate K-zz — so 
that counterpart of K.^z = A-z in the pseudo-Riemannian case is the pair of the equations 
K-^xx = ^2;y and K.yy = Ai-x- In particular, we have two real equations in both Riemannian 
and pseudo-Riemannian cases. 

Proposition 6.2. (a) After changes made by the rules described above, Theorems \l.l\ and 
li.i^ remain valid in the pseudo-Riemannian case. 

(b) Let g be a pseudo-Riemannian metric with vanishing ipi-, ^2 ^'^^ non-vanishing Then 
g admits a nontrivial linear integral. 

Proof, (a) The proof of Theorems ll.ll and ll.2l is done by means of formal calculus in which 
the condition of positivity of the metric was never used. 

(b) Considering in [Eil pp. 323-325] Riemannian metrics g on surfaces with non-constant 
curvature -R, Eisenhart introduces the following coordinate system: One of the coordinates 
is the curvature R itself and the other is orthogonal to the first one. In other words, g = 
gndx"^ + g2idy'^ and R{x, y) = x. The necessary and sufficient condition for the existence 
of such coordinates is that Vi? is not orthogonal to itself. In the Riemannian case this 
simply means the non-vanishing of Vi? and thus the non-constancy of R, in the pseudo- 
Riemannian case this is the condition ipi ^ 0. The rest of Eisenhart 's proof works also in 
pseudo-Riemannian. □ 
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6.4. Pseudo-Riemannian metrics with |V-R|^ = 0: local classification and nonexis- 
tence of cubic integrals. The goal of this paragraph is to prove 

Proposition 6.3. // a pseudo-Riemannian metric g with |V-R|g = admits a non-trivial 
cubic integral, then R is constant. 

We will need the following 



Lemma 6.3. Let g be a pseudo-Riemannian metrics with VR ^ and |V-R|g = 0. Then in 
generic null- coordinates it has the form 

f2{x)fi{y)y dxdy 
9 = ifi-,, ft 6.5 

with some functions fi{y), f{x), f^ix), and then the curvature is given by 

p _ 2/(x). 

Moreover, there exist null- coordinates in which g has the normal form 

9 = T^W^2 ^'ih R = 2f{x)^. (6.6) 

[y + f[x)Y 

where f{x) is some function. 

Such a metric g admits no Killing vector field, ie., no linear integral L. 

Remark. The generic formula is given up to exchange of null-coordinates {x,y). In (16.50 
and in the proof subscript indices fi{y)y,Uxyy, ■ ■ ■ denote usual (partial) and not covariant 
derivatives of functions fi{y), u{x, y), ■ ■ ■ 



Proof. Choose some null-coordinates {x,y). The condition |Vi?|^ = means that the gra- 
dient V-R is a null-vector. Hence R depends only on one of null-coordinates, say on x, 
and Ry = 0. On the other hand, R^ 7^ by the hypotheses of the lemma. Inverting the 
coordinate x if needed, we may assume that g = e^^^'''^^ dx dy for some function u{x,y). 
Then R = e~^Uxy and the equation |V-R|^ = reads [e~'^Uxy) = 0, or u^yy — u^yUy = 0. 
Substitution Uy = v yields Vxy — vvx = 0, which is equivalent to {2vy — v'^)x = 0. Hence 
2vy = f ^ + ip{y) for some function ip{y). This is a Riccati ODE, whose generic solution can 
be found as follows. Let vi{y) be a fixed special solution of the equation 2vy = v"^ -\- ip{y), 
so that ip{y) = 2{vi{y))y — fi(|/)^. Since ipi^y) was arbitrary, we may assume that vi{y) 
is an arbitrary function and ip{y) is given by the above formula. Now make the sub- 
stitution v{x,y) = vi{y) + Then the equation 2vy = + ip{y) transforms into 
Wy + Viw + 1 = 0. Since Vi{y) is an arbitrary function we may assume it has the form 
Vi{y) = '^j^^y"^ with an arbitrary function fi{y). Then the equation Wy + viw + 1 = reads 



fi{y)yWy + fi{y)yyW + fi{y)y = 0, which means {fi{y)yw{x,y) + fi{y))y = 0. Consequently, 
w{x,y) = ~^^^^f~^Y^ with an arbitrary function f{x), and 

n,\ _ h{y)yy _ o /i iy)y 

Now integration yields u{x, y) = log(/2(x)) + / f (x, y)dy which gives the desired form of the 
metric e'^^^''^^dxdy. 

The functions fi{y), f2{.x) can be easily eliminated using the substitution y' = fi{y), 
x' = Jf2{x)dx. 
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Now assume that a metric g = A(x, y)dxdy with Ry = admits a hnear integral L. This 
is equivalent to the existence of Killing vector field ^. Since the curvature R is constant on 
lines X = c, ^ must have the form ^ = a{x,y)dy. The Lie derivation gives 

=^g(A dx dy) = a dyX dxdy + X d^a dxdx + X dya dx dy. 

In particular, d^a must vanish. Then a depends only on y, and after the appropriate 
reparametrisation ^ = dy. But then A must be independent of x which would imply the 
vanishing of R. □ 

Proof of Proposition IKR We use the following notation: x, y are null-coordinates in which 
the metric has the form (16.61) . / = f{x), fi{x) = fi, f2{x) = f2 are functions of the variable 
X, f[,f2,f" their derivatives. Further, we set Ai = ai{x)dl, A2 = a2{y)dy, and then a[,a2 
denote the derivatives (in y in the case 02). The usual (non-covariant!) partial derivatives 
are denoted by d^^yK and so on. Besides we introduce the function Y{x,y) := y + f{x). 
Clearly, it satisfies the relations dxY = f and dyY = 1. Finally, ci, C2, . . . will be constants. 

In the case g = with Y = y+f{x) the only non-zero Christoffel symbols are = — ^ 
and = —y- This gives us the expansions 



2 2 ^ , 2 , Ya\ — Qf'ai 

,yy - JyyK + -dyK aud g^ = . 



K,yy = d'yyK + -dyK ^^1;^ 

Thus K.yy — g'^Ai .^^ appears to be an ODE 



Kix, y) = A(x) + ^ + (6.7) 



''^y + f {y + ff 

on K with respect to y with rational coefficients in which x is a parameter. The direct 
integration gives 

„,\ - f , Mx) , Ya 

and hence 

d.Kix,y) = fi + ^-g^ + 4,-i^-i^<±i^ + 'J^ (6.8) 
In the same way we obtain 

K,., = {d, + '-f)d,K and A2 = ^^^^ . 

Now we make the following observation: Substituting of (16. 7p in Y^i^d^^K + ^d^K^ 
and using dxY = f',dyY = 1 we obtain a polynomial (say P{Y)) of degree 5 in F whose 
coefficients are smooth functions in x, such that a2{y) satisfy the equation Ydya2 — Qa2 = 
P{Y). For any fixed x this is a polynomial ODE on 02. Integrating it, we obtain that (for 
any fixed x\) is given by 

a2iy) = ■ (/3(x) + / F-^P(F)rfr) . (6.9) 

In view of this formula, the problem of compatibility of the 3-codifferential A = {Al, A2) with 
the metric g = ^^_^^^'|^^^2 is equivalent to finding the functions f (x) , fi{x) , f2{x) , f3{x) , a{x) 
whose substitution in (16. 9p gives a function independent of x. 

Further, if we write P(Y) = X]^=o^'«^*' then the integration of (16.91) gives a polynomial 
expression for 02 

a2 = f,Y' - Elo 

with coefficients depending only on x. After expanding Y = y + f{x) we still obtain an 
expression for a2 which is a polynomial Q{y) = Yl^j=o V degree 6 with coefficients 

qj{x) depending only on x. However, since 02 is independent of y, all these coefficients must 
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be constant. In this way we obtain 6 equations qi{x) = Ci (ODE in general) on the above 
functions f (x) , fi{x) , f2{x) , fs^x) and ai{x). The expressions qi{x) can be calculated (or 
checked) using the formulas above, and we simply provide the final results. 

We are going to solve these equations using the prolongation (few times) and projec- 
tion(mostly) method, resolving subsequently equations qe{x) = C6,g5(x) = C5, . . . and sub- 
stituting the results into the next equation. The obtained new ODEs will be denoted by 

Eq, -^5, . . . 

The first equation is simply fsi^x) = Cq, so /a is a constant. The next equation g5(x) = C5 
(up to constant factor) reads as C5 + Gcg/ — /" = or 

f^ = c, + 6cj. (6.10) 

Substituting the latter expression q4{x) = C4 we obtain the next ODE E4 which reads 
C4 — ISce/^ — f[f' — — 5/C5. Resolving it we obtain 

= -2f[f - IO/C5 - 30c6/' + 2c4. (6.11) 

Next one substitution gives the equation which can be resolved as follows: 

<' = 6C3 + 120c6/=^ + 2/2/" + 60/2c5 - 24/C4. (6.12) 

Doing the same with q2{x) = C2, we obtain 

< = i2fy\- 2C2 - 3a[f - a J'" + I2/C3 + 6006^ + 40f C5 - 24/%). (6.13) 

The result of the same procedure for qi{x) = ci, qo{x) = cq is 

a[ = {Wf'Y' (30/'c3 + 60cef - UaJ"f' - AOfc^ + 50 f% + 5ci - IO/C2) , (6.14) 

ai = (2/'3)"i (2/ V - 2/='c3 + Co - 2fc, - 2c^f - f{x)c, + fc^) . (6.15) 

Recall that we have substituted in each successive equation the results of preceding cal- 
culations. This was clearly the projection procedure, and the next one is the prolongation 
one. Let us denote the expressions in (16.1211051) by ^3, . . . ,^0, so that the formulas read 
a[*'' = Ai. We compute solely the consistency equation Ai — d^Ao- The result is 

Ai - d,Ao = ji^{l6cef + 16/^C5 - 16/% + 16f C3 - 8pC2 + 8/ci - 8co.) (6.16) 

Since this consistency condition must be satisfied identically, we conclude that there are two 
possibilities: either /" vanishes identically, or all constants cq, . . . , ce must be zero. By (16. 6p . 
/" = means that the curvature R is constant. 

We exclude this possibility and consider the alternative case cq = . . . = ce = 0. Then from 
f l6.15p we see that ai vanishes identically. The equation fl6.12p reads now /2/" = 0, and so /2 
also vanishes. Finally, from (16. lip we obtain f[f' = 0, which means that fi{x) is constant. 

Summing up, we conclude from (16. 7p that the the function K{x,y) is constant. Finally, 
let us observe that the equation {y + f{x))dya2{y) = 602(2/) admit no non-trivial solution 
depending only on y. □ 
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7. Conclusion 

We presented an algorithm that, given a pair [g, A), where g is a. two-dimensional metric, 
and y4 is a holomorphic 3-codifferential, answers the question whether there exists a cubic 
integral whose Birkhoff-Kolokoltsov 3-codifferential coincides with A. Moreover, in the most 
interesting cases covered by Theorems II. ![ 11.21 it provides precise formulas for the integral. 
The algorithm works in an arbitrary coordinate system: we need to calculate certain ex- 
pressions given by precise algebraic formulas including the components of g and A and their 
derivatives, and compare them with zero. It is easy to implement the algorithm on modern 
computer algebra packages, say Maple® or Mathematica®. 

Our results suggest the following program for search for new natural Hamiltonian system 
admitting integrals polynomial in momenta integrals of degree ^ 3. 

Let us take a metric on a surface whose geodesic flow admits a cubic integral. For example, 
we can take a metric of constant curvature, or a metric admitting a Killing vector field, or a 
metric coming from one of the known natural Hamiltonian system via Maupertuis principle. 

Let K := Iplg he the kinetic energy corresponding to such metric and F3 be the integral. 
Let us now look for a function Fi : T*S' — >■ M linear in momenta, and for a function : S M 
such that F3 + Fi in an integral for K + V, i.e. , {F3 + Fi, K + V} = 0, where K is the kinetic 
energy corresponding to g. By Maupertuis principle, the later condition is equivalent to the 
statement that for every constant h the functions Hh := (^h-v) ^9 ^'^ ~^ ' 

commute: {Hh, Fh} = 0. 

Clearly, the Birkhoff-Kolokoltsov 3-codifferential of the pair Hh, Fh does not depend on the 
parameter h and coincides with the Birkhoff-Kolokoltsov 3-codifferential of the pair K, F^, 
i.e., is known since we know K and F3. Thus, the conditions Q2 = 0, Q3 = (or V = 0, 
^2 = O5 ^3 = case covered by Theorem II. 2p can be viewed as quasilinear PDEs 

on the coefficients of Fi and on V, i.e., on three unknown functions. Since the conditions 
are fulfilled for every h, the system of PDEs for V and Fi is overdetermined. Easy analysis 
shows that it is of finite type and there exists an algorithmic (though highly computational) 
way to find a solution. 

In other words, we suggest to look for new systems by adding potential to the known inte- 
grable geodesic flows, and the results of our paper ensures that one can do it algorithmically. 

Note that all known natural Hamiltonian systems either come from physics (e.g. Goryachev- 
Chaplygin), or were obtained by the naive version of the above procedure (for example those 
obtained in |Du-Mal IWTi 

As the most promising metrics g to start the above program, we consider the metrics 
admitting linear integrals. Recall that such metric are completely described. Local de- 
scription is due at least to |Dal §§590-593], global (= when the manifold is closed) can 
be found for example in [Bo-Ma- Fo] . Note that the starting point for the systems from 
|Du-Mal [Set IVe-Ts[ IYe[ were certain metrics admitting linear integrals. This will be the next 
direction of our research; we are quite positive that it is possible to describe all natural sys- 
tem admitting integrals polynomial in momenta of degree ^ 3 such that the corresponding 
metric admits a Killing vector field, and hope to find new examples. 

Another perspective application is to try to prove/disprove the following conjecture from 
[Bo-Ko-Fo] : if a real-analytic metric on the 2-torus admits a cubic integral, then it admits a 
linear integral, or its weaker form: if a real-analytic natural Hamiltonian system on the 2- 
torus admits an integral of degree three in momenta, then it admits a linear integral. Indeed, 
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as we explained in Corollary I2.1( ii). the space of holomorphic 3-codifferential on the torus 
is a 2- dimensional linear space. More precisely, there exists a global (periodic) coordinate 
system {x,y) on the torus such that the metrics has the form X{x,y){dx'^ + dy'^); in this 
coordinates system the form A is (c^ + ^i)^®^®^, where a, 6 G M. Then, one can view 
the equations ^2 = 0, ^3 = on the torus as two PDE on the coefficient A depending on two 
parameters a, b. 

One should also mention that in paragraph 12.41 we invented a calculus adapted to our 
problem. The calculations are much easier in this calculus (for example, in the recent papers 
[Kr] and [Br-Du-Ea] . which deal with a priori easier quadratic integrals, most calculations 
were done using computer algebra programs; in our paper, everything is done "by hand"). 
We expect that the same calculus will effectively work in the case of integrals of higher 
degree. 
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